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Abstract 



We propose, using the example of the 0(4) sigma model, a general method for solving integrable 
two dimensional relativistic sigma models in a finite size periodic box. Our starting point is the 
so-called Y-system, which is equivalent to the thermodynamic Bethe ansatz equations of Yang and 
Yang. It is derived from the Zamolodchikov scattering theory in the cross channel, for virtual 
particles along the non-compact direction of the space-time cylinder. The method is based on the 
integrable Hirota dynamics that follows from the Y-system. The outcome is a nonlinear integral 
equation for a single complex function, valid for an arbitrary quantum state and accompanied by the 
finite size analogue of Bethe equations. It is close in spirit to the Destri-deVega (DdV) equation. 
We present the numerical data for the energy of various states as a function of the size, and derive 
the general Liischer-type formulas for the finite size corrections. We also re-derive by our method 
the DdV equation for the SU (2) chiral Gross-Neveu model. 
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1 Introduction and Summary 



The study of the properties of Quantum Field Theories (QFT's) in finite volume, or at finite 
temperature, has a long history and numerous applications. Matsubara description [I] of finite 
temperature T thermodynamics, by considering the system in the periodic imaginary time t, has 
lead to the extensive study of the Euclidean QFT's with one compactified dimension with numerous 
physical applications [2]. 

Liischer found the leading finite size corrections to the mass gap in relativistic two dimensional 
QFT's [3, 4]. These corrections depend solely on the asymptotic S-matrix of the theory. Recently, 
Liischer corrections to various multi-particle states in integrable 2D QFT were conjectured [5]. 

For the integrable 2D QFT's, as understood during the last two decades, the ambitions can be 
much higher: these systems are usually solvable at any finite size though a systematic approach to 
such solutions, as well as a good understanding of the working prescriptions, are still missing. 




Figure 1: Physical channel, cross-channel and finite volume vs finite temperature. 

There are two main schemes to address the finite size calculations. The first, pioneered by 
Destry and deVega (DdV) [6], is based on the integrable discretization. Once such discretization is 
at hand, the system can be studied by the well established methods based on the transfer matrix 
approach and the resulting non-linear integral equation (NLIE), often called the DdV equation, 
calculates not only the ground state energy but also the spectrum of excited states. The method 
appeared to be very powerful when applied to the Sine-Gordon model [7, 9, 11, 12, 13], or to more 
general RSOS models [14], Toda theories [1-^)], hard hexagon models [16], etc. 

However, for generic integrable QFT it is far from easy to find the corresponding integrable 
lattice regularization and for many models such discretization is not known. Nevertheless, the 
problem can be usually tackled by using a computation scheme alternative to the DdV approach. 
As explained in the seminal work of Al.Zamolodchikov [17] this is achieved by the double Wick 
rotation trick: using the Matsubara imaginary time formulation we can first find the free energy in 
the infinite volume but finite temperature. Next we flip the meaning of euclidian time and space 
directions on the cylinder: r — > a, a ^ t, and interpret the free energy as the ground state of the 
system in finite volume L = (see fig.l). In this way we can obtain the exact finite volume ground 
state energy. This computational scheme is known by the name of Thermodynamic Bethe Ansatz 
(TBA). 

The TBA equations, whose number is infinite in many interesting models, can be usually con- 
cisely casted into the so called Y-system functional equations [18, 19]. Often the latter one can be 
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Figure 2: Plots of energies E of a few excited states of 0(4) model on a circle of a circumference 
L. The vertical axis corresponds to the values of -^E, the horizontal axis - to the values of L in the 
logarithmic scale. The lowest curve depicts the vacuum energy. The next one, labeled as Oq, shows 
the mass gap energy. The corresponding state is in the U{1) sector, with a single particle at rest, 
hence with the mode number = 0. The next states in the U{1) sector are denoted by 6nin2n3,-, 
according to the mode numbers ni,n2,n^, . . . excited for the 1-st, 2-nd, 3-rd, etc., particles. For 
all these states the SU{2)l and SU{2)r spins of the several particles are pointing in the same 
direction, say they are spin "up". The dashed line represents a state having a polarization out of 
the U{1) sector, with left and right "magnons" excited - it corresponds to the quantum state of 
two particles where both SU{2)l and SU{2)fi spins are in the singlet s = state. The qualitative 
explanation of these graphs will be given in subsection 5.2. 

rewritten in the form of DdV equations or some similar set of integral equations for a finite set of 
functions. The method was successfully used for many relativistic models [20, 21, 22, 23, 24]. As 
explained in the previous paragraph the computation of the exact ground state energy by means 
of this method is a relatively straightforward task with solid theoretical foundations. 

To obtain the exact spectrum comprising all excited states of the theory is, on the other hand, 
a much more involved - and a very interesting - task. A possibility to describe the excited states 
within the TBA approach, by modifying the analytical properties of the thermodynamic functions, 
was first suggested in [24]. Another possible way to obtain the spectrum of the theory, proposed 
around the same time, is based on the analytic continuation of the ground state energy with respect 
to the parameters of the model, such as the mass or the chemical potential, in order to find the 
excited states [25]. If the integrable lattice regularization is absent, it is not well understood why 
these methods work. Nevertheless, the results are usually in the excellent agreement with the 
perturbation theory, Liischer finite size corrections and the direct Monte-Carlo study for a wide 
range of sizes L (see for example [26, 27, 28, 29, 30] for 0{n) and related cr-models). 
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Figure 3: Domains of applicability of different descriptions of an integrable field theory at a finite 
volume L. In the ultra-violet regime, for small volume measured in units of a dynamically generated 
mass, the theory could be described by a conformal theory. In the infrared, at large volume, one 
can use the asymptotic Bethe equations. The leading order finite size corrections are governed by 
the (generalized) Liischer corrections. At any volume but for the ground state energy only one can 
use Thermodynamical Bethe ansatz. Hirota equation, equivalent to Y-system but more efficient 
when it comes to imposing appropriate analyticity properties, is a universal tool covering the whole 
diagram. 



For models with diagonal scattering, like the Sinh-Gordon theory [31], the whole classification 
of excited states is possible [32] . The situation is much more complex when we deal with the non- 
diagonal scattering. The nested structure of the corresponding Bethe ansatz equations leads to 
complicated magnon-type excitations and bound states. Little is known about the excited states 
in such finite size systems. The only models where the polarized excited states were investigated, 
using the DdV equations, are the Sine-Gordon model [7] and its supersymmetric version [11] as well 
as the tricritical Ising model [10]. By the existing methods only the sectors with diagonal scattering 
can be studied efficiently, as was done for example for the 0(4) sigma model in [30]. A general 
and unified description of all excited states of the cr-models like 0{n) or the SU (N) principal chiral 
field (PCF), and similar ones, having a "geometric" target space, is still absent. 

The main goal of the present paper is to give a method of a general and systematic description of 
all the states of integrable QFT's in finite volume. We will explain how to go beyond the asymptotic 
spectrum and compute the full finite size spectrum comprising all excited states of integrable sigma 
models. We do it here on the example of 0(4) sigma model and also for the SU{2) chiral Gross- 
Neveu model but our formalism is certainly more general and is most probably applicable to any 
integrable 1-1-1 dimensional cj-models. The main ingredients of the method are: 

• The two-particles S-matrix for integrable system allows us to write the periodicity condition 
quantizing the momenta of the physical particles on a large circle of length R. The equations 
following from the periodicity condition are so called asymptotic Bethe ansatz (ABA) equa- 
tions describing all states of the model. The details of this computation for the SO (4) sigma 
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model are given in Appendix A^. They are, however, vahd only in a sufficiently big volume 
compared to the typical interaction distance, Rm ^ 1 where m is the infinite volume mass 
gap. 

• For the ground state, the double Wick rotation (cr, r) — > (r, a) allows to reduce the problem 
to the thermodynamics. One can put the euclidian theory on the torus with one radius, 
R, very large and another one, L, arbitrary (see the fig.l). The ground state energy for a 
finite radius is related to the thermodynamic partition function. The exact equations for it 
can be found using the asymptotic spectrum given in the cross channel by the asymptotic 
Bethe equations. The resulting infinite series of integral equations, thermodynamic Bethe 
ansatz (TBA) equations, are casted into a functional form called Y-system. Here is the main 
assumption: we assume that different solutions of the Y-system describe not only the ground 
state but all the excited states. One should furthermore restrict the class of solutions by 
assuming certain analytic properties which will in particular identify the quantum numbers 
of the states we are considering. 

• Classical integrability of the Y-system, as a finite difference equation equivalent to the Hirota 
difference equation [,34, 19], allows us to express explicitly the infinite number of the unknown 
functions through a finite number of the basic ones [35, 36]. 

• The Baker- Akhiezer function of the Lax pair associated with the Hirota equation can be 
interpreted as the Baxter function encoding the "magnon" Bethe roots, responsible for the 
SU{2)ji and SU{2)l polarizations of states. The analyticity properties important for the 
full formulation of the resulting non-linear integral equation, are also suggested by the Lax 
equations. The gauge symmetry of Hirota equations allows to explicitly fix the final nonlinear 
integral equation (NLIE) for each state of the theory. 

The resulting equation can be studied in various limits (such as Liischer finite size corrections 
or small volume, conformal limit) or solved numerically in a rather efficient way. The fig. 2 shows 
some of our numerical results obtained from the new equation, plotting the energy of various states 
as functions of the volume. When the similar results are available in the literature the agrement is 
perfect. 

The general scheme elaborated in this paper on the example of the 0(4) sigma- model should 
be applicable to all integrable relativistically symmetric 2D QFT's. It should be also useful for 
the study of finite size effects when the system does not look explicitly relativistic but allows the 
S-matrix description and this S-matrix obeys the crossing symmetry, like the AdS/CFT S-matrix 
[37, 38]. Y-system and Hirota equations give a unified and powerful point of view at all this 
subject since they solve in an almost trivial way the "kinematic" part of the problem related to the 
representation theory, whatever is the symmetry or supersymmetry of the model [35, 39]. 

Our method based on Hirota equation, being exact for any finite size L of the system, reproduces 
well various limiting cases (see the fig.3). For the large L, the energies of the states are well described 
by the Liischer corrections [3, 4, 5]^. We derive them here for a general state with arbitrary 
polarization, which is also a new result, extending some hypothesis existing in the literature [5]. 
For small L, our results are well described by the theory of three free bosons, as will be discussed 

^Wc arc unaware of the existence of such derivation of the Y-system for the PCF in the literature 
^Actually, as we will see from our numerics, Liischer corrections work surprisingly well all the way until 
Lm ^ 1. 
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in the paper. The results for various low-lying levels, including the cases of non-diagonal scattering 
which are new, are summarized in the fig. 2. 

Our resulting NLIE can be brought sometimes to a form similar to the DdV equation. In the 
cases when the latter is available it can even coincide with DdV equation (an example of the chiral 
Gross- Neveu model is considered in our paper). It would be extremely interesting to understand the 
relation between the solution based on the integrable lattice discretization of [ io] and our proposed 
integral equations. Nevertheless we should stress that the real power of our method should be in its 
universality: it should work in all situations when the TBA equations in the form of the Y-system 
are available. 



2 TBA and Y-system for 0(4) sigma model, or SU(2) 
Principal Chiral Field 

The method we are proposing it quite general and we hope that a wide range of models could be 
solved using it. However for the sake of simplicity we will exemplify it on the SU{2) Principal 
Chiral Field (PCF), equivalent to the 0(4) sigma model. In section 4 we will also consider the 
SU{2) Chiral Gross-Neveu model. 



2.1 The Model 

The action of the PCF is given by the usual expression 



S^ = \ f dtdx {d^X,)\ iZ^XS = 1, (1) 

^0 J a=l 



whose target space is . It is equivalent to the SU{2) ® SU{2) principal chiral field (PCF) whose 
infinite volume solution was given in [41, 42, 43]. Indeed, by packing the fields Xi into an SU{2) 
group element h = X4 + i Yl^=i XjCTj with aj being the usual Pauli matrices, we can re-write the 



SpcF = -Tr2 I ^^^^ tT{h-^dahf . (2) 



action as^ 

The spectrum of this asymptotically free theory in the infinite volume consists of a single 

2" 

physical particle of mass m = Ae ^0 , where A is a cut-off. Its wave function transforms in the 
fundamental representation under each of the SU{2) subgroups. Al. and A.Zamolodchikov [33] 
proposed the exact elastic scattering matrix for such particles: 

5,.(^) = 5o(^)M^5o(^)|^, ^o(^)= vi + |)r(-j) ' 

where R{6) is the usual SU{2) R-matrix in the fundamental representation given by R{9) = 6 + iP 
and P is the permutation operator exchanging the spins of the particles being scattered. This 



•^In the AdS/CFT literature one usually uses a/A : 



47r 
'^0 
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S-matrix was established due to: (i) analyticity, (ii) unitarity, (iii) absence of bound states, (iv) 
crossing. In particular, (ii) and (iv) lead to the following identity 

5o(^ + ^/2)5o(^-^/2) = ^^ (4) 

on the scalar (dressing) factor. We can use this S-matrix to study the spectrum of N particles in 
a periodic space circle of a sufficiently big circumference L ^ m~^. The spectrum can be defined 
from the wave function periodicity condition 

n 5(^fc-^,)n'^(^'^-%)i^)=«""'''"''^"''^i*)' 
j=k+i j=i 

which quantizes the momenta of the physical particles. The asymptotic spectrum of the theory put 
on a large circle of length L is then given by 

N 

^ = ^mcosh(7r%) (6) 
i=i 

where 6j are solutions to the Bethe equation (see Appendix A for more details). In what follows we 
will measure all dimensional quantities in the units of m. Diagonalizing the periodicity condition 
(5) in the physical space by the usual methods (see an appendix in [39] for this model) we get the 
main Bethe equation 



e 



-iLsinhine,) T[ Q2rn n ^ Q«(^j + i/2) Qt,(% + ^2) 



k 

The magnon rapidities Uj and Vj are fixed by the auxiliary Bethe equations 
Qu{uj + i) _ (p{uj + i/2) Q^{vj +i) _ (p{vj + 



(8) 



Qu{uj - i) (l){uj - i/2) ' Qv{vj - i) 4>{vj - i/2) 
where 

Qw{x) = JJ(x - Wj) , for u; = n, v, (9) 

j 

and (t){x) = \{j{x - 9j). 



2.2 TBA and Y-system 

As we mentioned in the introduction, the ground state energy Eq{L) for arbitrary L can be com- 
puted starting from asymptotical Bethe ansatz in the cross-channel. For SU{2) principal chiral 
field this is described in detail in the Appendix A. The output is that the ground state energy is 
given by 

^o(^) = -lJ decosh{7T9)log{l + Yo) , (10) 

where Yq is one out of an infinite number of y-functions Yn with n G Z obeying the TBA-type 
equations 

logYn + Lcosh{7rx)6no = s*log{l + Yn+i)il + Yn^i), n = 0±l,±2,... (11) 
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with s = 2cosh(7rx) ^^^^ * denoting the convolution. If log(l^(3;)) for any n have no 

singularities inside the physical strip — 1/2 < Im x < 1/2 we can easily invert the operator s* to get 
simply = ea^^ + e~2^=: and these integral equations can be rewritten in a functional, y-system 
form 

y+y- = (1 + y„+i)(i + y„_i) , (12) 

supplemented with the asymptotic boundary conditions for large x 

y„^e-^"°'^''('^^)''"0 xconstn. (13) 

The superscripts it stand for shifts of the argument by ±i/2 ^, 

/± = /(x±i/2). (14) 

Eq.(12) has however many solutions and only one of them really leads to the ground state energy. 

SU,(2) U(l) SU,(2) 
-2 -10 17 

--o-o-»-o-o-- 

V e u 

Figure 4: Dynkin diagram (three central nodes) and its extension for the niagnon bound states 
(grey nodes) reEecting the structure of the Y-system. The central, black node corresponds to 
the U(l) sector excitations of the model (9-roots), the upper and lower nodes correspond to the 
more general states for magnon excitations for the SU{2)l wing (u-roots) and the SU{2)pi wing 
(v-roots). 

It is commonly believed that certain other solutions there describe the excited states [25, 44]. The 
energy of the A^-particle excited states is again given in terms of Yq but is modified 

1 r ^ 

E{L) = "2 / c?6'cosh(7r6l)log(l + yo) + ^mcosh(7r0j) , (15) 

where the extra terms are inspired by the analytic continuation in L and the points 6j [25] are 
singularities of the integrand in the first term 

yo(0,±V2) = -l, j = l,2,...,iV. (16) 

As we shell see, the last equation is nothing but the Bethe ansatz equation for physical rapidities 
modified at the finite volume. The last term in (15) is generated from the integral (10) by picking 
up the logarithmic poles (16). 

Our goal in this section is to make use of the integrability of the Y-system rewriting it in the form 
of classical integrable discrete Hirota dynamics. This allows us to write down explicitly a solution 
for all Yn in terms of a finite number of functions. Then we will restrict ourself to a certain sub-class 
of physically relevant solutions with particular analytic properties. The analyticity will allow us to 
fix the functions completely and parameterize all the physical solutions for the excited particle 
states in terms of a finite set of complex parameters, Bethe roots, restricted by supplementary 
Bethe equations reducing in the infinite volume to the usual Bethe equation. 



^We will often use even a more general notation, like + • • • + — y^^, _|_ or / ■ ■ ■ = f{x 

ik/2). 
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2.3 Hirota equations 

The y-system equations eq.(12) can be seen as a gauge invariant version of the so called Hirota 
equation or T-system 

Tk{x + i/2)n{x - i/2) - Tfc_i(x)Tfc+i(x) = $ ( x + i J ) ^ ( x - 4 ) • (17) 



It can be easily checked [19] that Hirota equation is equivalent to the y-system eq.(12) if we denote 

^ ^ $(x + «|)$(x-i|) ^ ' 

At first sight, this is just another trivial rewriting of the TBA equations, however the Hirota form 
appears to be particulary useful. Using Hirota equation we can also write 

, , rfc(x + i/2)Tfc(x - i/2) , , 

^ ' $(x + i|)$(x-i|) ^ ' 

Let us point out here an important fact. By evaluating the above equation for A; = at 0j it i/2 
where 6j is a zero of Tq we observe that 

To(^,) = ^ yo(%±i/2) = -l (20) 

which is the Bethe ansatz eq.(16). We will use this fact to associate zeroes of Tq with physical 
rapidities. 

Since Ifc(x) are real functions by their physical meaning (for ground state they are ratios of 
densities of complexes and of their holes, see Appendix A) we can restrict ourself to the case when 
Tk are real functions and $ and ^ are complex conjugated functions. 

Hirota equation (17) is integrable and has a Lax representation through the auxiliary problem 

[35] 

rfc+i(x)Q ( x + ij ) -Tfc (x - ^ ) Q (x + i^ + i ) = ( x + i J ) Q ( x - i J - 



Tk-i{x)Q[x-i^-i]-Tk[x--]Q[x-i-] = -^[x-i-]Q[x + i-] . (21) 



The compatibility of these two equations for the bi-vector of functions {Q(x),(5(x)} leads to the 
initial Hirota equation. Here Q is the complex conjugate function to Q. Note that if Tfc(x) are real 
functions then the second equation is simply the complex conjugate of the first one after shifting 
/c — > /c + 1 and x ^ x + i/2. Two particularly useful relations from this Lax representation are 

Ti{x) = To[x - 1/2) h$(x)- 



Q{x) ' ' Q{x) ' 

T_,M = ro(. + ./2)^-*wJ^, (22) 

Note that the first relation in (22) is a generalization of the famous Baxter equation usually written 
for the spin chains. We will see that in the infinite volume limit <&(x) = ro(x + i/2) and that 
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these equations reduce to the usual Baxter equation for spin chains, where Ti plays the role of the 
transfer matrix in fundamental representation for the magnons of the SUr{2) wing of the theory, 
whereas the second equation plays a similar role for the SUl{'2) wing (see Fig. 4). 

The main advantage of the Lax equations (21) is that they are linear in and we can easily 
express any in terms of Tq, ^ and Q in the explicit form 



Tk(x) 



+ Q ( x + i^^ ]Q(x 



(23) 



.k + l 



$ (x - + ij) 



^ Q (x - + ij) Q{x-i!^ + ij) 



This leads to a quite general and explicit solution of the y-system via eq.(18). A nice feature of this 
form is that one can efficiently analyze the L — > oo limit and reproduce the asymptotic spectrum 
described by BAE eqs.(7,8). This will be the goal of the next section. 

Hirota and Lax equations exhibit several important symmetries. First of all a discrete symmetry 
exchanging the n-wing and the u-wing (right and left SU{2)): ^ Y-k is induced by 



^ ^ , Q ^ Q , Q ^ Q 



++ 



(24) 



which will be quite useful for our further constructions. Moreover, both equations (17) and (21) 
are invariant under the gauge transformation 



Tk{x) - 




<I>(x) - 




l>(x) - 


9{x 


Q{x) - 





k 



k 



g X 



Tk{x) 



(25) 



To preserve the reality of Tk we should assume that g is the complex conjugated function to g. 
These transformations leave Yk{x) invariant. 

The general solution of Hirota equation (17) can be also presented in a determinant form [35] 



T^(x) = h{x + ik/2) 



Q{x + i^) R{x + i^ 



Q{x-i^) R{x-i^) 



(26) 



where h{x) is a periodic function: h^'^ = h{x + i) = h{x) and Q, R are two linearly independent 
solutions of the Lax equations (21) related by the Wronskian relation 



^{x) = h{x + i/2) 



R{x) Q{x) 
R{x + i) Q{x + i) 



(27) 



This determinant form will be very useful when we will formulate the general solution of the finite 
size PCF system for any state. It is not absolutely necessary to use it, but it simplifies some 
derivations. 
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2.4 Asymptotic Bethe Ansatz and Classification of the Solutions 



The main problem in computing the exact spectrum of the SU (2) PCF is to find the physical 
solutions to the y-system (12) or, alternatively, to the Hirota equation (17), i.e., obeying the right 
asymptotic properties (13). Their classification is a complicated task, especially when we want to 
take into account not only the excitations of U{1) sector but also the "magnon" type excitations 
of SU{2)i and SU{2)ji sectors. The goal of this section is thus to identify the large L solutions 
to the y-system (12). The discussion in this section is not completely rigorous since our only goal 
is to get an idea of how asymptotic Bethe ansatz (ABA) eqs.(7,8) appears from the y-system. 
Together with the expression (6) the ABA equations must appear from the large L asymptotic of 
exact solutions, as yielding the leading order value of the full spectrum. 

The main simplification in the large L limit is that Yq — > 0. From eq.(13) we see that Yq — > 
we are left with two decoupled chains of equations for > and k < [13]. 
For each wing we can introduce two sets of describing the corresponding solutions of the whole 
T-system: and such that yfc>o iYk<o )can be expressed in terms of (T^) by the formula 
(18). Then Yq = imphes 

rl\ = , r{' = . (28) 
Let us focus on since we can always use the wing exchange symmetry (24) to describe . 

We have to identify the solutions which will lead to the asymptotic spectrum described by the 
ABA. It turns out that in terms of Hirota potentials it is enough in this limit to make very 
simple assumptions, namely: 

• All TJ^^q{x) are polynomials at L ^ oo. We denote in this limit Tq(x) w 11 ~ ^j) = '/'(^)- 

• Q^{x) is a polynomial with real roots which we denote Q'^{x) = Y\j{x — uj). 

Then from eq.(22) we see that 

= r„"(. + i/2) and m.) = TS(. + ./2)Q''(.-iHTS(.-./2)Q-(. + .) 

From the polynomiality condition for T^{x) and T^{x) we get precisely the auxiliary Bethe equations 
eq.(8). 

Finally, we should note that eq.(7) for the physical rapidities Oj is also satisfied. This follows 
from imposing Y{){9j it i/2) = —1 for all zeros 9j of Tq*, see (20). At first sight, this seems to be 
impossible to satisfy since, as we noticed, Yq{x) is small. However this smallness appears because 
Yq is proportional to g^-Lcosh{-Kx) -^^j^igj^ indeed small inside the physical strip — 1/2 < Im x < 1/2 
but is of order 1 on the boundary of this strip. To impose this condition we must first compute Yq 
to the next order. 



From (12) at n = we get 



rpU+ rpV + rj-}U— rpV- 



^0 ^0 - )2 • (30) 

Defining S{x) = HjLi Sq{x — Oj) we have, from the crossing relation (4), S^~^ S = <\)l^'^^ , so that 

_ ( Tt{x)T-_^{x)S\x^il2) \^ ( T^{x)Tl,{x)S\x + t/2) y 
° ° " V <t>\x-^/2) ) \ 4>\x-^|2) ) ' ^'^^ 
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from which we can identify Iq up to a zero mode factor of i/q = e ^^"^^"^^ which obeys UqUq = 1- 
Such factor should be included into Yq to ensure the proper asymptotic (13). Thus we find 

Yo{x) e"^-^(-)rr(x)T-i(x) ^'|^ + ^^/^l . (32) 



Evaluating it at a; = 0^ — i/2 and using eq.(29) we get 

Qu{0k + i/2)Q,{e 
Qu{0k-i/2)Qv{ek-i/2) 



. ^, JLsinhM,) QujOk + i/2)Qv{ek + i/2) -p|- . X ,ooA 



which is nothing but the main ABA equation (7) for the middle node in fig. 4. We use here the 
notations 

Qv{x) = Q^{x-i) , Qu{x) = Qu{x) , (34) 

to make the u- and v-wings more symmetric. The advantage of these notations is that the wing 
exchange symmetry eq.(24) simply exchanges Qy and and in the large L limit they are real 
polynomials. 

Finally, since Yq[x) is exponentially suppressed for real x we can drop the integral contribution 
in (15) which leaves us with the energy as a sum of energies of individual particles, precisely as 
expected from (6). 

Notice that the Zamolodchikov asymptotic scattering theory is implicitly contained in the Y- 
system, as we see from the appearance of the scalar scattering factor 5^ in the formula (32). 



2.5 Probing the finite volume 



Now, having established the solution at infinite volume, we need an insight into the analytic prop- 
erties of T-functions in a finite, though large, volume. Let us find perturbatively the finite L 
corrections for the simplest vacuum solution which for large L corresponds to = = 1, (j) = 1. 
From eq.(23) one can see that for this case, to the leading order, ~ A; + 1 which implies for 
Yfe ~ \k^ -|-2|A;|. Thus we are looking for a solution in the form 



Yk = \k\'^ + 2\k\ + , fc = -oo,...,oo 



(35) 



where the first two terms in the r.h.s. are the trivial solution at L = oo, where as ~ Yq are 
small. We will see that the solution for the perturbation is unique under the assumption that when 
— > oo the perturbation goes to zero y/j — > 0. The linearized Y-system in the Fourier form is 



k 



k + 2 



s Vk+i 



k + 2 
Vk H ; — svk^i 



, A; > 



(36) 



where yk is the Fourier transform of yk and s 



2cosh(a;/2) 



is the Fourier transform of the kernel 



^ ~ 2cosii(7r6>) • yo = ^ is a fixed function. We see that this is a second order recurrence equation 
which in general has two linear independent solutions. Fortunately it can be solved explicitly.^ The 
general solution reads 



Vk 



k{k + l){k + 2) 



(k+2)\i 



k + 2 



Ci{uj) + 



e 2 



e 2 
k + 2 



■^One can use RSolve function in Mathematica to find the solution. 
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The needed solution satisfying % = Yq , = corresponds to Ci = Yq, C2 = 0. Making the 
inverse fourier transformation we get 

_ fc(/c + l)(fc + 2) f 1 

vr [Ax^ + k^ 4x2 + (A; + 2)2 * • ^"^^^ 

It can be easily checked that the approximate Tk yielding this solution through (18) are 

T;:_, = T^_,c.k + ^^^^*Yo , k>0. (38) 

and 

■^'-' = '+ 4(. + i/2''+V + l *^°- 

The iO in this expression can be dropped when computing Y/oo from (18) but is included in this 
expression so that (18) can also be used for k = 0, for more details see the discussion in the next 
subsection. 

An important feature of this asymptotic solutions for T^., which should persist at any L, is that 
it acquires two branch cuts at x E Ii±ik/2 when L 00/^ 



2.6 Exact solution for the vacuum 

We will now extend the solution found in the previous section to arbitrary L. First, we notice that 
the solution in terms of Tk is much simpler than in terms of Y^. For the vacuum we can use the 
following ansatz inspired by eq.(38) 

Tk-i = k + ^^^^*f, A: = +0,l,2,... (40) 

where / is some function which for large L becomes Yq. One can easily see from the linear system 
eq.(21) at Q = Q = 1 that this ansatz solves the Hirota equation and can be presented in the form 
eq.(23) with ^(x) = Tq{x + i/2 + iO). Thus the Y-system equations eq.(12) for \k\ > 2 are satisfied 
automatically. Notice that none of the T^'s has singularities on the real axis, which is of course a 
necessary feature of the solution: the physical quantities Y^ should not be singular there. 

To check that the equation for A; = 1 is also satisfied we have to define Yq in terms of T^. For 
that we can simply analytically continue eq.(40) to the point k = +0 which gives T-i{x) = f{x)/2. 
We also have $(2;) = TQ{x-\-i/2 + iO), ^(x) = Tq{x — i/2 — iO) as mentioned above. These properties 
are supported by the second equation (22) which can be viewed as yielding the spectral density in 
terms of a jump on any of two infinite cuts. Then we get 

Y (^) = To{x + i/2-iO)To{x-i/2 + iO) _ ^ ^ Ti{x)f{x)/2 

°^ ' ro(x + i/2 + iO)ro(x-i/2-iO) Tq{x + i/2 + iQ)TQ{x - i/2 - ^ ' 

This equation relates Yq and /. With Yq so defined the Y-system equations at = 1 are now also 
satisfied. However the equation (11) for A; = is still not used. Using 



(i + yi)(i + y_i) = (i + yi)' 



Ti{x + i/2)Ti{x - i/2) 
Tq{x + i)To{x - i) 



^The term "branch cut" is not very appropriate here since the infinite cut has no branch points. However, 
as we shall see, a spectral representation will allow us to define Tk{x) in the whole complex plane in terms 
of spectral density integrals along the cuts. 
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and recalling that s is the inverse shift operator we obtain^ 

Mx) = e-^-^(- \ J^f (42) 
[To[x + i)To[x 

Combining it with eq.(41) we get 



. To (x + i/2 + iO)To(x-z/2-iO) eosh(..) 



which, in virtue of the eq.(40), gives a closed equation for f{x). 

Notice that from eq.(43) r_i(x) = f{x)/2 is exponentially small for large L with 

r_i(x) ~2e~^'=°"^('^^). (44) 

The finite L solution to equation (43) can be easily found by iterations, starting from this large 
L asymptotic and gradually diminishing L. We solved this equation numerically and get a perfect 
match with the existing results (see the Tab.l comparing our results with [26]). 



L 




Leading order 


Eq.(43) 


Results of [26] 


L = 


4 


-0.015513 


-0.015625736 


-0.01562574(1) 


L = 


2 


-0.153121 


-0.162028968 


-0.16202897(1) 


L = 


1 


-0.555502 


-0.64377457 


-0.6437746(1) 


L = 


1/2 


-1.364756 


-1.74046938 


-1.7404694(2) 


L = 


1/10 


-7.494391 


-11.2733646 


-11.273364(1) 



Table 1: We solve numerically eq.(43) the use Yq from eq.(41) to compute the energy of the 
ground state using eq.(lO). 

In the next subsection, we generalize this solution to the excited states in the U{1) sector. 



2.7 Generalization to U{1) sector 

In this section we will study in detail the U{\) sector of the theory where we consider the states 
with N particles with the same polarization, i.e. with no magnon excitations. Hence we can put 
all Q = 1. As mentioned before - see eq.(20) - for N particle states we expect T^iOj) = for each 
of N rapidities of the particles 6i, . . . ,9^ ■ 

In the previous section the vacuum state, with no particles excited, was analyzed. We saw that 
To{x) inside the physical strip, <I>(x) above the strip and ^>(x) below the strip could be described 
by a single function 

^(-) = 1 + 4^ * T-i , (45) 

such that 

{4>(x-i/2) ,Im(x)>l/2 
To{x) ,|Im(x)|<l/2 . (46) 

^{x + i/2) ,Im(x)<-l/2 
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Figure 5: The function F{x) in (47) can be recast as a contour integral as in (49) with the 
contours as represented in this figure. 



Here we build a generalization of (45) for the case when Tq has an arbitrary number of zeroes 
inside the physical strip for which (46) holds: 



F{x) = (pix) 1 



1 



1 



1 



(t>{y-^/'^)x-y + i/2 (t){y + i/2) X - y 



I ^ T^im\ (47) 



with (t){x) = njLi(2; — ^j)- The overall factor of 4){x) appears because T(){dj) = 0. The spectral 
representation of F{x) as two integrals over the two infinite cuts at Im (x) = ±1/2 is inspired by 
(45) and can be also seen from the linear problem (21). Indeed, we have 



r_i(x) = To{x + i/2) - $(x) = To{x - i/2) - l>(x) 



(48) 



which justifies the choice of spectral densities used in (47). To see that (46) indeed holds we write 
(47) as 

F(x) = Mx) f / T^^y)/^iy) I / dy ^y-i/2)/<l,{y) ^ f dy ^y + i/2)/m \ .^gs 
\ jty 27ri y — x jty+ 2TTi y — x J^- 27ri y — x J 

The contours 7, 7"^ and 7" encircle respectively the physical strip, the region above the strip and 
the region below the strip, see figure 5. For this relation to be equivalent to (47) we require that for 
large x we should have Tq[x), (^[x — i/2),^[x + i/2) — > (j){x) at \x\ — > 00 along the corresponding 
contour. Finally, for (46) to hold, the ratios in (49) should be analytic inside the corresponding 
contours. Notice that at large L the function r_i is exponentially small and thus Tq, ^>~'' — > </>(x) 
as expected from our discussion in section 2.4, to get the ABA equations. The large x limit should 
be similar to the large L limit since the source term in the y-system Q-Lcosh{-Kx) jg gj-^g^ii both 
cases. 



^Wc introduce a natural notation g* 
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Let us now consider the other Hirota functions T^. From (21) we have Ti(x) = Tq{x + i/2) + 
= <I>(x) + Tq{x — i/2) which in terms of the function F{x) reads 

Ti{x) = F{x + i/2 + iO) +F{x-i/2 + iO) =F{x + i/2- iO) + F{x-i/2- iO) , (50) 

so it is indeed regular on the real axis. Notice that Ti is regular at least inside the enlarged strip 
|lm (x)! < 1. In the same way we can easily see that Tfc>o is analytic inside the strip |Im {x)\ < 

Having expressed Tq, <I> and Ti in terms of r_i through the function F{x) we can find a closed 
equation of T_i from the Y-system equation for n = 0. The derivation is parallel to the one in the 
previous section and it leads to 

r_,(.) = (F(x + i/2) + F(. - i/2)) ^'-;'f+;°'f-:'g.-'°> , 

[F[x + i)F[x — i)\ 

supplemented by the quantization condition lo(^j + ^2) = —1. As before, the solution to these 
equations can be easily found from iterations as is explained in the Sec. 5. The numerically calculated 
energies of a few states of this U{\) sector are presented on the fig. 2. 

In the next section, we generalize these results to any excited states including the magnon 
polarizations. We will use a different strategy and incorporate the gauge invariance of Y-system to 
find the solutions of Y-system eq.(12) matching the L = oo asymptotic of the Sec. 2. 4. 



3 Finite size spectrum for a general state of PCF 

We will now describe how to construct the solution for the most general state of the PCF at finite 
volume L, having an arbitrary number of physical particles with arbitrary polarizations in the 
SU{2)ii and SU{2)l wings (characterized by left and right "magnons" Ui and Vi). Our method is 
based on the following observations and steps: 

• We know from eq.(32) the structure of the poles and zeroes of all Y^'s in the limit L — > co 
when Yq = 0. We assume that this structure will qualitatively persist even for finite L, and 
the classification of the appropriate solutions of the Y-system will follow the same pattern of 
poles and zeroes. 

• We will recast the Y-system in terms of T-system (Hirota equation) since the analytic struc- 
ture of Tfc's is much simpler than of as we saw from the vacuum solution (38) at L — > oo. 

• For any "good" solution of Y-system there is a family of solutions of Hirota equations related 
by gauge transformations (25). Hirota equation can be solved explicitly in terms of Tq, ^> and 
Q as in eq.(23). 

• For L — > oo we have two independent solutions for T^'s as we saw in the previous section. 
For one solution are asymptotically polynomials for /c > and for another one with 
/c < are polynomials when L is large. We can then smoothly continue these two solutions 
to finite -L's using the gauge freedom to preserve polynomiality of Q's. 

• We have two global solutions of Hirota equation which can be parameterized by T^J*, 

, and by Tg , Q^' . They represent however the same and unique solution of the Y-system 
and thus should be related by a gauge transformation g : Tg = g o Tf, see (25). 
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• Using certain assumptions about analyticity of Tq and supported by the Lax equations 
(21), we can express them as different analytic branches of the same analytic function Gu- 
The same can be done for Tq and —^v m terms of G^. 

• The solution will be completely fixed by the existence of such a gauge transformation g{x) 
which relates its u- and the f-representations. At the end we will have one single non-linear 
integral equation (NLIE) on g{x). 



The final equation for g{x) is new for the Principal Chiral Field. It is different from the system 
of 3 DdV type equations used for the same model in [30]. Still it resembles in many aspects the 
non- linear Destri-de Vega (DdV) equation which appears when studying other integrable models. 
Indeed, our method is very general and it allows to generate DdV- like equations for large classes of 
sigma models in a systematic way. For the models for which a DdV equation is known we expect 
our integral equation to coincide with it after an appropriate change of variables. We check this 
hypothesis on the SU (2) chiral Gross-Neveu model for which we re-derive indeed the known integral 
equation. 



3.1 Exact equations for the finite volume spectrum 

In this section, we will derive the finite volume spectral equations of the previous section in the 
most general form, valid for all excited states of the model with any number of physical particles 
with arbitrary polarizations (i.e. with any quantum numbers). 

As we discussed below in the infinite volume, the solution of Y-system with Yq = can be 
described in terms of two independent sets of Hirota potentials and . Since these two different 
solutions of Hirota equation correspond to the same solution of Y-system they are related by 
a gauge transformation g{x). These two solutions of Hirota equation can be continuously and 
unambiguously deformed all the way from very large L, where we know the solution (see the 
previous section), to any finite value of L. The gauge ambiguity for any of the two solutions, TJ^ or 
T^, can be fixed by choosing Qu and to be polynomials for any L. Of course we can no longer 
assume and T^, as well as the corresponding and to be polynomials. Instead we will 
assume certain analytic properties for them and we will see their consistency with the solution we 
find at the end. 

We introduce a polynomial (j){x) with real zeroes 9j, j = 1,2, . . . , N of Tq. They correspond to 
the rapidities of physical particles on the circle. The gauge function g{x) relating the two solutions 
of the T-system is assumed to be regular and to have no zeros on the physical strip, so that 
Tq = ggTQ has the same zeroes as Tq there. We also assume that*^ 

* <j>^x+i/2) ( <p[x-t/2) ) regular for Im x > —1/2 (Im x < 1/2) in the whole upper (lower) half 
plane and goes to 1 for |x| — > oo in all directions in the upper (lower) half plane; 

T^(x) 1 1 

• -^^^ is regular and goes to 1 at x — > ±00 inside the physical strip — ^ < Im x < i; 

^When a; ^ 00 we know that ^0(2;) — *■ ^ i.e. it is exponentially small, as in the case of 

large L, and the Y-system, as well as the T-system splits in to two independent u- and ti-wings with 
To(x) - $(a;) - 0(a;) - . 
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The first property is somewhat similar to the forth property from the previous subsection: the 
large x asymptotics is governed by the same exponential g"-^™*^'^'^^ as the large L asymptotics. As 
a consequence of the second assumption, inspired by the integral representation (47) for the U{\) 
sector, T^^q(x) are regular for — (/c + l)/2 < Im x < {k + l)/2 . 

Similarly, for another solution we assume that 

* 0(a+^72) ( cf,(x-f/2) ) regular for Im x > — 1/2 (Im x < 1/2) in the all upper (lower) half plane 
and goes to 1 when — > oo in all directions in the upper (lower) half plane 



is regular and goes to 1 at x 



iboo inside — i < Im x < i strip. 



Note that with these properties all Ifc/o are automatically analytic in the physical strip, as we 
see from (18). Yq is, in the strict sense, only analytic on the real axis, but the detailed analysis of 
the Appendix D shows that we can expect its analyticity even in a finite strip around the real axis. 
As concerns the T- functions, although we use in different circumstances TJ^ or T^, we will get the 
same result for all Y^, since they are related by a gauge transformation g{x). But analyticity will 
be explicit only for one wing at a time: at A; > 1 for TJ^ and at /c < 1 for T^. 

The listed properties are enough to explicitly relate , with the corresponding and 
using a certain integral representation for them. In Appendix B we follow this line of arguments to 
formulate the complete set of equations for an arbitrary state, including a NLIE for g{x) and the 
associated equations for the Bethe roots. However, it appears to be more advantageous, especially 
for the numerics, to use the integral representations for the logarithms of T^, T^, <I>^, We will 
derive in what follows the corresponding equations defining the energy of a general state. 
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Figure 6: The ratios indicated close to each contour are analytic inside the corresponding contour. 
Thus we can obtain them in each of these regions using a single resolvent Gu{x) as in (52). 



Let us define two new functions analytic on a Riemann surface with two infinite cuts at Im x 
±i/2 

2m y_oo x-y-i/2 2m x-y + i/2 



Gui{x) 



dy, {w = u,v) 



(52) 
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with the following spectral densities along the infinite cuts^ 

'T^{x + i/2) 



Pu{x) = log 

and their complex conjugates 

Pu{x) = log 

Then from (18) we have 



p^{x) = log 



Pv{x) = log 



-^v{x) 



-<^v{x) 



log(l + Yq) = Pu + Pu = Pv + Pv ■ 



(53) 



(54) 



(55) 



When L is large enough we know from the results of the section 2.4 that Tq{x) ~ ^wi^ ~ V^) — 
i>^(x + i/2) ~ (j){x) and thus to the leading order pw{x)^s are exponentially small. It is also 
important for our analyticity assumptions listed above to hold that pw{x) ~ g-^-coshTra; ^ ^ 
along the real axis. Together with these analyticity assumptions the following formulae are true at 
any L 



G,{x) 



log 
log 
log 



-^^ix-i/2) 
-^^ix + i/2) 



Imx > +1/2 
|Ima;| < 1/2 
Imx < -1/2 



Guix) 



log 
log 



log 



^u{x + i/2) 



Imx > +1/2 
|Imx| < 1/2 
Imx < -1/2 



(56) 



These formulas are easily understood from simple contour manipulation as depicted in figure 
6. Let us consider the resolvent Gu, plug (53) and (54) into (52) and consider separately the 

^"^(a:/^^ ^^'^ "^"^la)^^^ . Since — > 1, x ^ ±oo, we can close the 



terms containing 



4>{x) 

'^u{x-i/2) ^g^£|-gj. ^YiQ obvious shift of integration variable) around the upper 



contour in the integrals containing around the physical strip and contracting it around 

the pole y = X we obtain the middle relations in (56) if x lies in the physical strip. Similarly, 
using the fact that in the upper half-plane "^"[^^"'^^^ — > 1, x — > oo, we can close the contour in the 

integrals containing log — 
half-plane and contracting it around the pole y = x we obtain the upper relations in (56) provided 
Im(x) > 1/2. The integrals containing log ^^^^^^y^ are treated similarly with the contours being 
closed in the lower half-plane. For the resolvent G^ the same sort of reasonings apply. 

As we mentioned above, the two solutions of Hirota equation we defined in this way, are related 
by a gauge transformation = g o T^. However, the polynomials Qu and Qy are not necessarily 
related by this gauge transformation g{x). Instead, one can easily see that is mapped to another 
linearly independent solution of eq.(21) = We can use eqs. (26,27) to express all T^'s and 
<^'s in terms of Q's. In particular, we have 



h+ 



9'^ 



9^ 



rpU 



h 



0+0+ o^o- 



9 



9 



(57) 



^These spectral densities denoted by p should not be confused with the densities of Bethe roots g used in 
appendix A in the derivation of the Y-system ground state equations. 
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Similar relations for v wing can be obtained from the gauge transformation = g g^^u a.nd 
Tq = ggT^. For the densities (53-54) this yields 

Note that one can get from pu by exchanging indices u <^ v and 5 ^ 1/^. 

We see that the densities and thus both G„ and Gy now can be expressed solely in terms of 
three polynomials Q„, 0„, (f) and a function generating the gauge transformation relating the 
two wings. It is left only to find a closed equation on g{x). We do it in the following subsection. 

3.1.1 Closed equation on the gauge function g{x) 

In the previous subsection, we managed to express all relevant quantities in terms of three poly- 
nomials g^, gu, (/) s-iid a complex function g{x). Using the condition that two solutions of Hirota 
equation are related by the gauge transformation generated by g{x) we can write a closed equation 
on that function. In particular, using the fact that and are related by the gauge transfor- 
mations (25) we obtain 

$„ = (/V^n. (59) 

It gives a closed relation on g which we can rewrite, assuming that g is regular within the physical 
strip, as follows 

g ^ i^^iLsinhinx) , (60) 

where the zero mode of the inverted operator was chosen to ensure the proper large x asymptotic. 
Indeed with this choice T^^ = — -|t ~ g-icosh{7ra:) ig^ding to the right behavior of Yq (see eq.(18)) 
at large L. Using (52) this can be re-casted as 

g = ie^^^"'°^('^^)5(2;) exp [s * G^{x - i/2 - iO) - s * Gu{x + i/2 + iO)] (61) 

where we used (56) and the identity 

f_y' = (^S+S-y' = S , (62) 



following from the crossing relation (4). We remind that (j){x) = Y\f=i(x — 6j) and S{ 



X 



The closed NLIE (61) for g{x) is our main result. Together with the expressions for the densities 
in terms of g (53-54) it allows us to calculate (1 -|-lo) and thus to obtain the energy of a state (15). 

In what follows in this section we will rewrite it in a little more convenient form which will be 
useful for the numerical computations for particular states. 

Conjugating the last equation we find 

g = -ze-*^/2sinh(^x)^-i(^) *G^{x + i/2 + iO) - s * Gu{x - i/2 - iO)] . (63) 

Finally it is useful to translate these equations into an equation for the phase g/g 
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where Kq = log5g and we used 

Gw{x + i/2 + iO) - Gy,{x - i/2 - iO) = Ki{x + i/2 - iO) * - Ki{x - i/2 + iO) * 



(65) 



and the convolution form of the dressing kernel as Kq = 2s* Ki, Ki(x) = . For the squared 

norm gg we get from eq.(60) 



99 



n-iU~{~ rpU- 

^0 ^0 



rj-\V-\~rj-iV- 



rpv 



exp(G„ - Gu) 



(66) 



where we used y^Yq = 1 inside the square brackets to get the last equality. This equation can 
be also obtained from the gauge transformation Tq = ggT^ . As we shell see in section 5 these 
equations can be efficiently solved numerically, by iteration, where at each iteration step a single 
convolution integral arises involving the densities pu and p^. 

We use eq.(58) and eq.(59) together with our analyticity assumptions to constrain g,^u,v and 
Q ^. In the next section we will fix the remaining finite number of complex parameters - zeros 
of polynomials Qv, Qu and the real zeroes of <j). After that one can use eq.(23) to construct all 
Tfc. In appendix C we show that all obtained in this way will be real functions and thus Hirota 
equation for them is satisfied. This means that we solved indeed the Y-system with the right 
physical analytic properties for the solutions. 



rpU,V 



3.1.2 Finite size Bethe equations and the energy 

Finally it is left to explain how to fix the finite number of constants, the Bethe roots 9j, Uj and Vj, 
which are zeros of the polynomials (j) , Qu, Qv and which completely characterize a state. The zeros 
of <f> are by definition the zeros of Tq which means that at these points e''"^^^^*/^) = e''"'-^^^*/^-* = 
as we can see from eq.(58). It can be also written as follows 

±}L±lLl = 1 at x = 0., (67) 

which we can rewrite using eq.(64) as 



e-^'^^'-'^^-'^'^ = -52(^,)|^ exp Q [K^ * (p^ + p.) - K+ * (p„ + p.,)] ) , at x = 0,- . (68) 



Note that when L — > oo we can neglect the last factor to get precisely the usual infinite volume 
ABA eq.(7). The equations for the auxiliary Bethe roots Uj can be derived in many alternative 
ways. The most standard way is to demand analyticity of Ti at x = Uj (see eq.(22)) 

^uQu'^ + To"-Q++ = , at X = Uj- . (69) 

We see that in general there is no reason to assume Uj to be real when L is finite. Using the resolvent 
Gu to represent Tq and appearing in this expression we get the auxiliary Bethe equations 
following from (69) in the form 

1 = P(n.) , (70) 

where P{x) is defined on the upper half plane by 

P{x) = e-x.^[Ki{x-i/2)* Pu- Ki{x + i/2)*pu] , Im x > , (71) 
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In the large L limit P{x) ~ 1 and we get the ABA eq.(8). A similar equation fixes the roots Vj. 

The integral equation (61) together with the equations (68), (70) fixing the zeros of the polyno- 
mials Qu, Qi,, 4> are the complete set of equations which one should solve to find the full spectrum 
of the SU{2) X SU{2) principal chiral field. Once g{x) and the positions of the zeros 9j,Uj,Vj are 
found, we can compute the exact energy of the corresponding quantum state from eqs. (15,55) and 
(58) 

N ^ 

E = cosh(7r6'fc) " 2 / cosh(7r2;)(pu + Pu)dx . (72) 



k=l 



where we can use due to (58) Pu + Pu = Pv + Pv = log 
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Let us remark that our construction for a general state in this paper was based on the assumption 
that in the asymptotic regime L — > 00 all the roots Uj^Vj become real. However, it is well known 
that the complex solutions are also possible. We hope that even in this case our equations maintain 
their form, although this situation deserves a special care. 

In section 4 we will explain how to efficiently implement these equations for numerical study. 
Before that, in the next section we will study the large L behavior of these equations thus repro- 
ducing not only the large volume results of section 2.4 but also the subleading corrections (Liischer 
corrections) . 



3.2 Large volume limit: ABA and Liischer corrections 

The SU (2) principal chiral field spectrum is given by (15), or (72). As we have seen in the previous 
section, in the large L limit the Bethe roots 6j are given by their asymptotic values obtained from 
a solution to the asymptotic Bethe equations, and since Yq and p's are exponentially small we can 
drop the integral contribution in (68) and (70) and recover the usual asymptotic spectrum. In this 
section we focus on the leading finite size corrections to this result. 

Due to these corrections auxiliary roots Uj and Vj become complex even if they were real 
asymptotically at large L. In this section we denote the real part of the roots Uj,Vj by Uj and Vj 
while the (small) imaginary parts we denote by Atij and Auj. The positions of the momentum 
carrying roots 6j are also corrected, however they stay real. We will also use the notation 

^u{x) = n - ' ^''(^) = n - • (^3) 

i 3 

We have to compute the first correction to the positions of the Bethe roots. To the leading order 
we can drop the exponentially small densities p^ and p^ in eq.(61), to get 

(7(x) ~i5(x)e^^/2sinh(^:.) (74) 

which we can use to compute the spectral densities from (53-54). We see that some terms in the 
expression for pu are exponentially suppressed and we can expand 

^ Qu-Qu , -Lcosh(.x)'g+ Qr + Qu-Qu , ^t^T^l 
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In the last equality we neglect the small imaginary part of the axillary roots and we use (22) to 
the leading order together with the gauge transformation T'_^i = g~g~^T'_^i and the crossing relation 

s+s- = C. 



The poles at x = Uj should cancel, due to eq.(70), among the first and the second term since 
the density by definition is regular. We introduce the notations pl^^ and pu^ for the first and the 
second term in (75). The first one can be simply written as 



(76) 



^ X — Uj 



Since the whole density is regular we can apply the principal part prescription to the finite integrals 
in (64) without changing the result. Having done so we are free to split the convolutions into 
convolutions with p^^^ and p^'^\ In (68) we should then expand the factor 



^ exp (^Ira j {x - y)p^^^ (y)^ exp (^Im j {x - y)pf^ (y) 



(77) 



and the similar factor for the v roots, to the next to leading order. We notice that pu^ is purely 
imaginary to the leading order, as seen from eq.(76), and therefore we can simplify the term in the 
square brackets 

Si „„„ ( 1 , , . m\ ^ « /, ^ . „(i)\ _ ci A , rf'* + p« \. at 



where the convolutions are understood in the sense of principal value. Thus in the Bethe equations 
(67) in this approximation the imaginary parts of the axillary roots cancel against the contribution 
from p(^) and we simply get 



QuQv ^ 



at X = 9j . (79) 



Proceeding in the same fashion in the eq.(70) for the auxiliary Bethe roots we arrive at a similar 
conclusion. Namely only the real parts of the auxiliary roots survive when we separate the density 
into p^^^ and p^^^ 

/ „-r (2) 



= exp ( -2ilm A7 * p\f>] at x = Uj . (80) 

See appendix E for details. We see that all terms except for the convolutions with p^'^^ have a simple 
effect of absorbing the imaginary parts of the Bethe roots. 

It turns out that the remaining convolutions, containing p^'^\ can be nicely written in terms of 
the leading order Yq found before in (32), 



Qu J \ Qv 

where the 6j appear in and S while the uj (vj) auxiliary roots appear in the corresponding Baxter 
polynomials Qu (Qv)- Notice that this quantity is already exponentially small, so we can take here 
the asymptotic values for the auxiliary roots. As explained in detail in appendix E, the quantities 
inside the principal part integrals are related to the derivative of this function with respect to 6k or 
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Uk and Vk which we treat in (81) as independent variables. More precisely we have the remarkable 
identities 



^lrn{K-ie.-y)[pi^\y)+£\y)]) = (82) 

mm[K^iu,-y)pi^\y)) = +^^^ • (83) 
Thus we finally obtain the corrected Bethe ansatz equations in the following elegant form: 



.^.Lsinh(..)^2|^|^ ^ «^P(/^^^^?^) atx = 0,, (84) 



— xjt — exp -/" — ^ c^V at X = 

(j)- £1++ ^ \ / 9^,- ^ / J 

It is not completely surprising that we managed to express everything in terms of Yq. To the leading 
order, Yq can be expressed in terms of S-matrix only: it is the relevant eigenvalue of the operator 

g-Lco.h(..)^^ (5oi(x - 0i)5o2(x -92)... Son{x " 9^)) . (85) 

We see that (84) corresponds precisely to the conjectured equation (27) in [o] only inside the U{1) 
sector. However, our result is different from outside the U{1) sector when there are axillary roots 
Uj and Vj. Finally the equation for the energy of the state corrected by the finite size effects is 
given by eq.(15) in terms of Yq in the leading approximation and the corrected positions of the 
roots 9j which should be found from eq.(84). 

The right-hand sides of the corrected Bethe equations (84) have a simple interpretation: for the 
middle equation, it reflects the contribution of scattering of the "physical" particles off the virtual 
ones on the cylinder, whether as the other two reflect the same effect for the "magnons" responsible 
for the isotopic degrees of freedom of the particles. Although these equations are derived here only 
for a particular model their form looks very universal and can be immediately generalized to any 
other integrable sigma model where the exact scattering matrix is known. 

3.2.1 Single particle case 

In this section we consider the single particle case for the Liischer-type correction of the previous 
subsections. This analysis was done in a more general context in [5]. 

When we have a single particle with momentum 9i (84) yields simply 

Lsinh(7r0i) = 27m - j ^de,YQ{y) , (86) 

which corrects the leading order quantization condition 

Lsinh(7r6'?) = 27rn. (87) 
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Now, from (81) we see that the x dependence in Yq{x) comes from the exponential factor e~-^'^™^'^('^^) 
and also from the combinations x — 9j,x — Uj and x — Vj appearing in the remaining terms in this 
expression. Thus 

Ju Jv 

dyYoiy) = -L7rsinh(7ry)yo(y) -^de,Yo{y) -J^duMy) -J^^vMy) , (88) 

k=l k=l k=l 

which, in the case we are considering, with = 1 and = Jv = 0, allows us to simplify (86) to 

Lsinh(7r6'i) = 27rn + J yLsinh(7ry)yo(y) , (89) 
so that the leading finite size correction to the energy (15) reads 

E{L) - cosh(7r05') - j cosh(7ry) (l - tanh(7ry) tanh(7r05')) e-^™'^i>(^s^)tr Soi{y - 0?) , (90) 
precisely as expected for the Liischer corrections [4]. 



4 SU (2) Chiral Gross-Neveu model and related models 

Our NLIE resembles the Destri-deVega equation and, at least in the cases the last one is known, can 
even coincide with it. In the cases when the DdV equation is not known, like the SU{2)l x S'C/(2)/j 
PCF, or 0(4) model studied in this paper, we obtain a new, DdV-like equation. In this subsection, 
to demonstrate our method, we show how to reproduce the DdV equation for the chiral SU{2) 
Gross-Neveu model on a finite circle. 

The TBA equations for this model are given by the same Y-system (12) with an important 
difference that Ys<o = (see Fig. 8). In particular, since y_i = we have T_2 = which implies 
due to the eqs.(21,17) 

= Trf - ^Q~^ , T+{rZi = ^-^+ . (91) 
Then it is immediate to check that the quantity 



satisfies, due to eq.(91), the relation 

Note that ^ is a pure phase on the real axe. Thus, restoring the proper zero mode exponential, we 
find 

n+ /<t>\ 

As before, to make this a closed equation on A we introduce the resolvent G and the density p 



27ri J x — y — i/2 27ri J x — y + i/2 ' $ 
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Again analytic properties of Tq and $ lead to eq.(56). Using the linear problem (21) we can write 
r_i(5++ - T+Q = and we see that 



log 



Q 

Q 



(96) 



Thus A satisfies the closed equation for A since <^ = G{x + i/2 + iO). We see that A{x) plays a 
similar role as g{x) in PCF. We can easily compute Iq in terms of A 



l + Yo = ePeP = (1 + A+) (l + IM") 



(97) 



We see that the factors of Q cancel from this expression. It is also possible to write eq.(94) in a 
simpler form without factors of Q using the useful identities. First using eq.(65) and assuming that 
the density is regular not only on the real axes but also for — 1/2 < Im x < we write 



where p~+ = p{x — i/2 + iO). Using the following relation 



P 







= log 


Q 


[q- 




stQ-. 



-i^o * log 



where 

Seix) = So{x - 0j) , Suix) = So{x - uj) , Suix) = So{ 
j=i j=i j=i 

we can get rid of Q's in (98) and finally obtain the known DdV equation 



N 



X — Ui 



A 



(98) 



(99) 



(100) 



(101) 



Then for Yq we get the standard relation which should be used to compute the exact spectrum from 

1 



E{L) 



-— / m 
2 



cosh(7r0) log(l + Yq) + m cosh(7r^j 



(102) 



Furthermore, both BAE - for the physical rapidities 9j and for the magnon rapidities uj - can be 



written 



A{dj) = -1 , A+{uj) = -1 . 
We also notice that A{x) has poles at x = Uj + i/2. To see that we use 

ro(0,) = o ^ 1 + ^(0^.) = S^S^ = o. 



$(nj)~l ^ l + ^+(uj) 



_ Q{uj)T+{u,) 
Q{uj) ^Uj) 



0. 



(103) 

(104) 
(105) 



Equations (100), (97), (102) and (103) are precisely the DdV equations derived in [6, 7]^°! As shown 
in this section our method can be directly generalized to other models whose TBA y-system 



^°In these papers the sine-Gordon model was considered. The SU{2) Chiral Gross-Neveu model is a simple 
limiting case of this theory, see discussion in the paragraph below. 
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equations are known. It would be very interesting to make a systematic study of such models using 
our formalism. 

For example, a simple generalization of the case studied in this section is obtained by considering 
the functions T^, $ as well as the Baxter functions Q to be periodic in the imaginary directions 
with period iv. This amounts to considering the trigonometric solutions of Hirota equation (17) 
and the corresponding linear problem (21) - this should correspond to the sine-Gordon model [13]. 
We take the Baxter polynomials Q{x) and the large L limit of Tq, <I>^, ^ ~ (/> to be 

^ -J-r smh-ix — Uj) ~, , SL- smh-ix — dj) 

Q(^) = n -^^^ ' ^(-) = n 



instead of polynomials. Then most of the formulae in this section go through with minor modifi- 

' x-i/2 ^ 
x+i/2 



( X ill \ 

cations. For example, instead of the SU{2) Chiral Gross-Neveu S-matrix So = ( ^A/o ) we will 



find the sine-Gordon dressing factor 

i 



smh|^x-i/2) _ . psin(u;x) sinh(^L^) 
'-'0 — I . , TT / — ; — ^TTvT ~ — zexp 



sinh^(x + i/2)7 w cosh(|)sinh(fw) ■ 

Thus it seems that our method allows to derive the sine-Gordon DdV equations of [6, 7] in an easy 
way. For ^ oo we recover the SU{2) chiral Gross-Neveu model. For an integer v the Y-system 
can be truncated as represented in figure 8, see e.g. [13]. 

Another interesting class of models which one could analyze using our formalism is represented 
by the so called sausage model (see e.g. [45]). This model can be considered as a generalization of 
the 0(4) model, or SU{2) PCF, in the same sense as the sine-Gordon model is a generalization of 
the SU{2) chiral Gross-Neveu. The inhomogeneous XXX-spin chain present in the SU(2) Gross- 
Neveu model and describing the isotopic degrees of freedom is generalized in sine-Gordon model to 
the XXZ chain, with the anisotropy parameter u introduced above. Similarly, the sausage models 
can be seen as two interacting inhomogeneous XXZ chains parameterized by the inhomogeneities 
6i,. . . ,9n and anisotropies i' and u' . It would be very interesting to generalize our 0(4) model 
results to this more general class of models. 

Our approach to deriving DdV like equations is strongly based on a smooth interpolation 
starting from the IR asymptotic Bethe ansatz description; hence, by construction, our states are 
very well identified. On the other hand we did not carry out a detailed study of complex solutions 
such as the states represented by Bethe strings in the large L limit; for these states some of 
our formulae might need to be modified. Within the DdV approach based on descritizations of 
integrable models, many interesting complex solutions were studied: e.g. holes, special objects, 
wide roots, self-conjugate roots etc. It would be interesting to complete our approach to include 
all physical complex solutions and thus obtain a precise dictionary between these two approaches. 
In particular this would teach us which solutions to DdV NLIE correspond to physical states. 



5 Numerics 

In this section we explain how to efficiently solve numerically the equations derived in the Sec. 3. 
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5.1 Implementation of numerics and Mathematica code 



For simplicity let us focus on the U{1) sector where we consider M-particle quantum states with 
M spins pointing in the same direction in SU{2)l and SU{2)^. Then our equations simplify 
considerably as was explained in subsection 2.7. First of all, since there are no spins excited we 
have trivial Baxter polynomials Q„ = Q„ = 1. Thus, from (58) we see that Pu = Pv = P with 



log 



1 



1 



(107) 



We also notice that since this is a symmetric configuration where the u and v root configurations 
are the same (there are no roots at all) we have, see e.g. (66), gg = 1 and thus g{x) is a pure phase. 
In particular, for real x, we can simplify the density to 



p{x) = log 



1 



1 



from where we see that we can express it solely in terms of g~^ 
only to determine its argument from (64) which now reads 



(108) 

Since 5 is a pure phase we need 



Kn *l0g 



1 



\s 



+ |2 



1 



(109) 



This is almost perfect for numerical implementation but still needs to be slightly improved. The 
reason is that we want to iterate this equation by evaluating the right hand side for real x. But this 
will yield the updated values of g{x) in the left hand side whereas for the next iteration we would 
need g~^{x). To fix it, we simply shift x ^ x + i/2 in this equation and define A{x) = {g~^{x))^ to 
get 



M 



A 



-L cosh(7ra;) 



J]5o2(x-0,+i/2)exp Ko*log 



A-l 



\A\ 



1 



K++ * log 



A-l 
\A\ - 1 



log' 



A-l 

A\ - 1 



(110) 

where the convolution of Kq~^ is understood in the principal part sense. We have explicitly written 
S{x) to render the presence of the Bethe rapidities more explicit. These are fixed by the main 
Bethe equation (67) which in our notations is simply 



M 



X 



6j ) exp I 2i Im 



Kq *log- 



A-l 



1 , X 



:iii) 



For completeness let us present here the Mathematica code to solve these two equations by itera- 
tions^^ It is a slightly simplified, and thus less efficient, version of the code we used for the plots 
in figure 2. 

First we introduce the S-matrix and the kernel Ko, 



S0[x_]=I*Gamina[-(x/(2I))]Gamma[l/2+x/(2I)]/(Gainma[x/(2I)]Gamma[l/2-x/(2I)]) ; 
KO [x_]=D[Log[S0[x]-2] ,x]/(2*Pi*I) ; 

^^One can copy the code directly to Mathematica from .pdf 
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Next we specify the size L and the mode numbers n = {rei, . . . ,njvf}. For example, if we want 
to study the system with L = 1/2 and three particles with zero mode numbers we write 

n = {0, 0, 0>; 
M = Length [n] ; 
L = 1/2; 

We will perform several integrals from — oo to +cxo but the integrands have exponential tails so 
that it is quite useful to introduce a cut-off X for all the integration intervals at this point. A 
reasonable cut-off is given by ^-LcoshiiTX) _ iq-8 Furthermore, at each iteration step we will have 
to construct an updated function A{x) which we do by means of an interpolation function, 

X=ArcCosli [8Log [10] /L] /\ [Pi] ; 

F [S_] : =Funct ionlnterpolat ion [S , {x , -X , X} , InterpolationPoints->30] ; 
Next we introduce 

eq[i_ , v_] : =L*Sinh [Pi*x [i] ] +Sum[If [i==j , ,Log[SO [x [i] -x [j] ] '2] /I] , {j ,M}] -2n [ [i] ] Pi+v [ [i] ] ; 
BAE[v_] :=Table[x[i] , {i ,M>] / .FindRoot [Table [Re [eq[j ,v] ] ,{j,M}] , 
Table [{x[i] ,2i/M-l/2},{i,M}]] ; 

which yields the solution to Bethe equations Lsinh(7r0j) -|- "^Zj^i J ^^S ^oi'^i ~ ^j) ~ '^'^iT-i + Vi where 
Vi is a perturbation to the equation number i. Comparing with (111) we see that this perturbation 
at step k will be given by the convolution appearing in (111) evaluated at the solution 9i computed 
in the previous step. The leading order BAE's correspond to = and are thus given by 

\ [Theta] [0] =BAE [Table [0 , { j , M}] ] 
Also to leading order the function A(x) will be simply given by 

A [0] =F [-Exp [-L*Cosh [Pi*x] ] *Product [SO [x-\ [Theta] [0] [ [j] ] +1/2] ~2 , {j , M}] ] 

Then we introduce the density p and its conjugate p at the A:-th iteration step as 

r[k_,y_] :=Log[(A[k] [y] -1) / (Abs [A [k] [y]]-l)] ; 
r c [k_ , y_] : =Conj ugate [r [k , y] ] ; 

Finally the code 

A[k_] : =A[k]=F [-Exp [-L*Cosh[Pi*x]] Product [SO [x-\ [Theta] [k-1] [[j]]+l/2] ~2,{j ,M>] 
Exp [Nlntegrate [KO [x-y] r [k-1 ,y] -KO [x-y+1] *rc [k-1 ,y] +1 , {y , -X, x ,X} , 

Method->Principal Value] -2X-rc [k-1 , x] ] ] ; 
phase [k_] [x_] : =Nlntegrate [21m [KO [x-y-1/2] r [k-1 , y] ] +1 , {y , -X , X}] -2X ; 
\ [Theta] [k_] :=\ [Theta] [k] =BAE [Table [phase [k] [\ [Theta] [k-1] [[j]]] ,{j,M}]] ; 

yields the k-th iteration quantities in terms of those computed at the {k — l)-th step. The energy 
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Figure 7: Difference between the predictions from tiie asymptotic Bethe ansatz and from tlie 
generalized Liischer formulas to the exact (numerical) results for the two particle state polarized in 
both SU{2) [u,v,9-ii]. The solid (blue) line represents ^<'^'^^~^aba ^Ii^Iq ^jjg dashed (red) curve 
depicts ^exact-EL^scher _ jg clcar that the latter approximates the exact results with excellent 
accuracy, especially for large L of course. 

of the state is then given by 

En[k_] :=Sum[Cosh[Pi\[Theta] [k] [ [j] ] ] , {j , M}] -NIntegrate [Re [r [k,y] ] Cosh [Pi*y] ,{y,-X,X}] ; 
For example, to obtain the result of the first 8 iterations we simply run 

Table [En [k] , {k, 0, 8>] 

to get {10.2414,10.2425,10.2425,10.2424,10.2424,10.2424,10.2424,10.2424} where we notice that the 
iterations are clearly converging to the exact value 10.2424 up to the precision we are working at. 
It is instructive to compare these results to the value predicted by the asymptotic Bethe equations 
alone, 

Ebae = 10.3388 

and to that predicted by the generalized Luscher formulae discussed in section 3.2 which gives 

-E'Luschcr = 10.2396 

5.2 Discussion of numerical results 

Now we will try to interpret the behavior of various states on the fig.l as functions of the volume L. 
Let us start from the vacuum, the lowest plot there. At very small L, the 0(4) model should become 
a 2d CFT of three massless bosons: if we introduce in (112) a rescaled field Y = ^(L) (Xi, X2, X^) 
and X4 = Y^l - e^{L){Y^^ + + Y^), where eg(L) ~ is the effective charge, very small in 

this limit (the effective radius of the sphere R{L) = ^(L) is very big), the action will be 
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S^= dtdx ^{dM^ + O (e^(L)) . (112) 

a=l 



In the ground state, the Casimir effect wih define the limiting energy: Eq ~ — 1^+0(1/ log^(L ^)), 
with the central charge c = 3, which gives Eq^' — "I? the value compatible with —0.18 of the 
fig.l. 



The energies of excited states are 

N 



k=l a=l 



where = (ra^^\ n^^\ n[.^^) are the momentum numbers of particles constituting the state. We see 
that the small L asymptotics of our plots are well described by this formula: The excited states 
in the U{1) sector, denoted in the fig.2 by 9ni,n2n3,-, according to the mode numbers ni, 712,713, 
approach the values predicted by (113) (up to the circumstance described in the last footnote). In 
this sector they have no excited left and right magnons (no u,v roots), and only one component is 

activated: j^Enmins-i^) = "-i^^- Say, the curves 6^0, 6*00, 6*000, 6*0000, approach -1/4 at L 0, 
the curves 9i, 9qi, ^qoi approach 3/4, the curve 62 approaches 7/4, etc. The state with one left and 
one right magnon excited, denoted as [u,v,9^i^i], also approaches 7/4. 

The qualitative behavior of the states 0o, ^00, ^000, ^0000, etc, at very small L's can be explained 
by the fact that the quantum fields are dominated by their zero modes. ^'^ Since the momentum 
modes are not excited the field Y{a, r) does not depend on a. The action and the hamiltonian 
become: 



e^(L) J 4 L 

where the angle (P{t) represents the coordinate of a material point (a top) on the main circle of the 
unit sphere, and J is the corresponding angular momentum. The quantum mechanical spectrum 
of this system is well known: 

L , „ „ , 1 o , , , , 771(771 + 2) , , 

m times 

where a; = 7 + log (^'^ ^^ti^^ ^ — 13.66 is a constant [-[('>]. This formula explains well the fact that 
the corresponding plots on fig.2 converge slowly, as inverse logarithm, to —1/4 and their spacing is 
approximately linearly growing with the number 777,. 

The perturbative calculation of the mass gap [EQg{L) — Eq{L)] for L <^ 1 was done in [40] and 
was compared with the numerical results following from the TBA approach in [26] . Since our 
numerics is in a perfect agreement with [26], for the states for which their method works, we will 
not review it here. We only recall that, in the logarithmic approximation. 



^^Thc convergence to the limiting value is very slow at L ^ 0. At L ~ 0.1 for our calculations we are still 
far from the limiting value of the energy. In ['2()], where the numerics reached L = 10^^, the result is —0.226, 
considerably closer to the limiting value. 

^■^We would like to thank A.Tsvelik, P.Wiegmann and K.Zarembo for the explanations on this subject. 
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^[Ee.{L)-E,{L)]^\^^, (L«l) (116) 

which is in the perfect agreement with eq.(115) at m=l. We also compared eq.(115) for m = 2 
with our numerics and found a good agreement^^. The inverse logarithm in eq.(115) explains well 
the slight divergence of various curves with zero mode numbers at increasing L in fig. 2. Most 
probably, the divergence of the other plots at increasing L, corresponding to the same value of 
Y^k=i Sa=i I'^i"'*! be also perturbatively described by the same inverse log terms. It would be 
interesting to study the small L limit analytically to recover analytic properties of the perturbation 
theory. 

At large L, we enter the realm of the asymptotic Bethe ansatz with Liischer-type exponentially 
small corrections to these regime. Actually, they describe very well our exact numerical plots for all 
the considered states, considerably beyond the values of L allowed by the approximation, as seen 
in the plot 7. 

6 Conclusions 

We derived in this paper the non-linear equations yielding the energy of an arbitrary excited state 
in the 0(4) two-dimensional sigma-model, equivalent to the SU{2) principal chiral field, defined 
on a space circle of an arbitrary length L (measured in infinite volume mass gap units). The main 
formulae we found are reviewed in subsection 6.1. 

Although we considered mostly the 0(4) sigma model the new method which we develop here 
is very universal and should be applicable to any integrable relativistic sigma model, such as the 
SU{N) principle chiral field at any N , 0{n) sigma models at any n, or more exotic models like 
SS-model or supersymmetric Sine-Gordon model (see [29] for the examples and the ground state 
energy and [11] for some excited states). In figure 8 the ground state Y-system diagrams for many 
known models are represented ~ it would be extremely interesting to perform a systematic study 
of such models using our formalism. 

We also hope that the method will eventually allow to calculate the spectrum of finite size 
operators (such as the Konishi operator) in the N=4 SYM theory, when applied to its dual, the 
integrable string sigma model on AdS5 x background, on a world sheet cylinder of a finite 
circumference, as inspired by the works [38, 47, 48, 49]. In spite of the last spectacular applications 
of the S-matrix approach for the perturbative calculation of wrapping interactions for Konishi and 
other twist-2 operators [5], the problem of finding the dimensions of such operators at any coupling 
is still open. 

The Hirota equation, which is equivalent to Y-system, appears to be a remarkable tool for 
solving the integrable sigma models in 2 space time dimensions. Not only does it help to collapse 
an infinite system of equations into a few ones, but it also helps to guess the analytic properties of 
the remaining unknown quantities and thus formulate the problem in terms of a single equation for 
a complex function g{x) . This function has the transparent meaning of a gauge transformation 
between the T functions in two different solutions to the Hirota equation, but solving the same 
Y-system. This equation reminds of the famous Destri-deVega equation, and in the known cases, 

i-^The discrepancy with the r.h.s. of eq.(115) for m = 2 and L = 1/10,1/100, 1/1000 is 0.057,0.034,0.023. 
^■"^It might be a few functions for other sigma models but always a finite number of them. 
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like the chiral Gross-Neveu model, even coincides with it in certain variables, as we demonstrated in 
this paper. However, for many interesting sigma models the Destri-deVega equation is not known 

- in particular for general states in finite volume systems. Our method suggests a systematic way 
of deriving such DdV-like equations. For example, in the case of SU{N) symmetry, we can expect 
that the closed set of such equations should not contain more than 4 real functions (or two complex) 

- the total number of the gauge functions for the general Hirota equations. It would be interesting 
to apply our method to the SU (N) principal chiral field [41, 42, 43], especially in the large limit, 
which is explicitly solvable for a non-zero magnetic field [50]. 

There could be interesting applications of our method to conformal QFT's in two dimensions 
if we consider them as some limiting cases of massive theories. These limiting cases could be the 
ultraviolet limit of a small volume L ^ oo (see for example [14]) or analytic continuation w.r.t. the 
number of components of a field, like in [51], or something else. 

Finally, one of the most promising grounds for the applications of our method should be the 
case of supersymmetric sigma models, a quickly developing subject, which is very useful in many 
physical problems ranging from AdS/CFT correspondence to disordered systems. The method of 
solution of Hirota equations applied for the supersymmetric spin chains with the symmetry algebras 
gl{K\M), was worked out in [52, 53, 39, 54]. 



6.1 Summary of the main formulae 

In this subsection we summarize our final integral equations in a self-consistent set of formulae. 
The main example considered in this paper was the 0(4) sigma model where the particles have two 
SU{2) spins as internal degrees of freedom. To compute the exact energy of A^-particle states with 
Ju left spins down (and thus N — left spins up) and right spins down (and thus N — Jy right 
spins up) we should solve the single integral equation on a complex function g{x) 



g[x) = iei^^^'^^^''^S{x) exp [s * G„(x - i/2 - iO) - s * Gu{x + i/2 + iO)] , 

.■ r(i/2- 

'r{l/2+ix/2)T{-ix/2) 



where * stands for convolution, s{x) = 2 cosh -kx ' '^(^) ~ Hj^i Soix — Oj), Sq{x) = ^x/2)^(+^x/2) 



and the resolvents are given by 

where the densities are parameterized in terms of g{x) as 

gPii — .2 fsPv — y, 

'^Qt^Qv~ -QuQv ' '^Qt^Qu- -QAJ 

with Q«,(3;) = nfc=i(^ ~ ^fc) and Q«,(x) = nfc=i(^ ~ '"^fc) being the Baxter functions encoding the 
Bethe roots of the left and right "magnons" {w = u^v). The superscripts it indicate the shifts by 
±z/2, so that e.g. Q^^ = Q^ix -|- i) and the bars indicate the complex conjugation. Finally the 
constants Uj and Vj are fixed by the finite volume Bethe equations 

Q+(^^)Q+(^5(%) = 1 r(^,)Q++(%) r{vi)Qt^{vi) , 

Qu{ej)Qy{e,)gi9j) ' 0+(n,)Q--(n,-) ' <A+(^j)Q^"(^j) ' ' 



34 



where Pw=u,v{x) is defined on the upper half plane by 

2 1 

Pw{x) = exp [Ki{x - i/2) * p^, - Ki{x + i/2) * p^] , Ki{x) = — 



vr 4x2 + 1 ' 

and by its analytic continuation in the full complex plane. The energy of the state is then given by 

E = cosh(7r6'fc) - / cosh{-Kx){pu{x) + pu{x))dx . 
k=i •'-'^ 
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A Derivation of the (ground state) y-system 

To compute the ground state energy Eq{L) of the SU{2) principal chiral field in a periodic box of a 
size L we can compute its Euclidean path integral Z with the fields living on the space-time torus 
of the size L x R , where the periodic imaginary "time" R is very big 

Z = e-^^o(^) . (117) 

Following Al.Zamolodchikov [17] we can compute this quantity exchanging the role of L and R so 
that the latter becomes the space variable whereas the former becomes the periodic time. Since 
i? — > oo the spectrum corresponding to the new Hamiltonian can be computed from the asymptotic 
Bethe ansatz and the finite periodic time L means that we should consider the system at a finite 
inverse temperature L. Thus we conclude that 

Eo{L) = f{L) . (118) 

where f{L) is the free energy per unit length of the SU(2) PCF at the temperature 1/L in the 
thermodynamical limit, when R ^ oo. 

To compute the free energy we will start by reviewing the asymptotic spectrum of the theory 
as given by the asymptotic Bethe ansatz. Then we will recall what are the magnon bound states 
(complexes, or strings) and how they are organized in the complex plane. We will see that the 
quantum states in the thermodynamic limit can be described by the densities of these complexes 
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and their holes. From this description we will write the entropy formula and thus find the desired 
exact free energy as the result of the saddle point approximation at i? — > co. This will give the 
TBA equations. 

Particles in the SU{2) principal chiral field transform in the bi-fundamental representation under 
two SU{2) groups. The theory is integrable and thus the general scattering process factorizes into a 
sequence of two-body scattering events. The 5-matrix [33] describing the scattering of two particles 
with momenta and energies given by 

Pj = msmh(TT9j) , Ej =mcosh(7r0j) , (119) 

depends only on the difference of rapidities 9 = 9i — 62 

Su{9) = 5o(^)M ^ So{9)l^^ , M9) = , (120) 

where R{9) is the usual SU{2) R-matrix in the fundamental representation given by 

R{9) = 9 + iP, (121) 

where P is the permutation operator exchanging the spins of the scattered particles. 

From now on, we will measure the length L in the units of the mass gap m, which means that 
we will put m = 1. 

When particles are put on a large circle of length R the periodicity condition to be imposed 
on the wave function reads 

-f{9j) e*-R«i>^h(^e,)^ ^ ^ ^ ^^22) 

where T is the transfer matrix 

f{9) = tro (^01 {9-9^)... Son {9 - 9n))) , (123) 



with the index for an additional auxiliary particle which we scatter against all physical particles. 
The trace is taken over this auxiliary space. Indeed, when the transfer matrix is evaluated at a value 
of the physical rapidity 9j the corresponding S'-matrix Soj{9 — 9j) becomes simply —Poj^Poj which 
means that the auxiliary particle changes the quantum numbers and becomes the physical particle 
9j. Then (122) becomes the periodicity condition (5) which physically states that once we pick the 
particle j and carry it around the circle the total phase acquired by the wave function - which will 
be given by the free propagation Rpj plus the phase shifts do to the (factorized) scattering with 
each of the other particles - must be a trivial multiple of 27r. 

Using the algebraic or analytic Bethe ansatz technologies it is possible to diagonalize T{9) for 
any value of 9 using the same eigenvector basis (see an Appendix from [55] for the details). Multi- 
particle states with left spins down (and thus N — left spins up) and right spins down (and 
thus N — right spins up) are parameterized by J„ auxiliary Bethe roots Uj and roots Vj and 

'^(^)^ = ^^^^"(^ - - V2)^ , (124) 

where T^^^'^ is the transfer matrix in the fundamental representation associated with the left (right) 
SU{2) spins, 

^ Qn,i9 + m9 - i/2) + Q^i9 - ^9 + i/2) ^ ^^^^^ 
Qw{9) 
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and 



N N Ju Jv 

m^UiO-Oj) , Sie) = llSoi9-9,), QuiO) = llie-Uj), Q,{9) = l[{9-vj). (126) 
j=i j=i j=i j=i 

The rapidities 9j and Uj,Vj are then fixed by a set of nested Bethe equations. The Bethe equations 
for the physical rapidities 9j are given by the periodicity condition (122) which can be written as 

g-imiisinhM,) ^ __^^M_rpu^0. _ i/2)T^{9j - i/2) , (127) 



or simply 



Qu{9,+i/2)Q,{9j+i/2) 



The magnon rapidities Uj and Vj are fixed by the auxiliary Bethe equations 

Qujuj + i) _ (pjuj + i/2) Qv{vj+i) _ (p{vj+i/2) 
Qu{uj - i) (p{uj - i/2) ' Qv{vj - i) (j){vj - i/2) 



(129) 



which appear in the diagonalization of the left and right transfer matrices. Notice that these 
equations ensure that the apparent poles in (125) drop out and render Ti{9) polynomial as it 
ought to be. For each solution to these equations we obtain the energy of the corresponding 
quantum state from 

TV 

E = ^cosh{TT9j) . (130) 
i=i 

To be able to compute the free energy f{L) we need to understand how the solutions to these 
Bethe equations organize themselves so that we can introduce the entropy density. Let us consider 
the auxiliary roots u, obviously the same considerations will apply for the v roots. These roots can 
take complex values. When uj has a positive imaginary part the r.h.s of the Bethe equations in 
(129) diverges, 

(piUj + i/2) N~*oo 



4>{uj — i/2) 

which means that 

Qujuj + i) 
Qu{uj - i) 



oo , 131 



(132) 



must diverge as well. This can be achieved if there is another magnon rapidity Uk such that 
Uj — Uk — i- Thus, in the thermodynamical limit the magnon rapidities will organize themselves 
into a Bethe-string of n roots Uj spaced by i. In particular, a single real root corresponds to a 
Bethe string with n = 1. The Bethe equations can then be multiplied for u^'s belonging to the 
same string so that this gives new Bethe equations, solely for the (real) center of each string. This 
is the usual fusion procedure applied at the level of the Bethe equations. The resulting equations 
look as follows. Introduce the magnon bound states: 

(n) ('^) I • ^ / I 1 \ ■ 1 

j,a = + i-[n + 1) - la, a = 1, . . . , n. 
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Multiplying the equations for a given n-bound state we get for eq.(129) and eq.(128) 

Pj^a j,n "a Uj 1 2 

(n) r. . -n M (m) -n+m (") ("^) -In— ml "+J^ , ^ ^ 

vy^9p + ij^ ^ u)' -ul ' - u)' -ul ' - 1'—^ 2 / u)' -ul ' + IS 

-'■-'■„(")_ _,■ 11 ,,("■)_ _L n+m „,W „,M I An-m\ li i „,(")_ „,M „ 

In the thermodynamic limit we will have a large number of each type of Bethe roots which we 
can describe by a density Qn- We use n = for the density of 9 particles, n > 1 to describe the 
density of u Bethe strings of size n and n < —1 for the v Bethe strings made out of — n roots. 
For each density of particles we also have the corresponding density of holes Qn- Bethe equations 
in the thermodynamic limit, obtained by taking the logarithmic derivatives of both sides of these 
equations, read 

R °° 

gn + Qn = — COsh(7r6l)(5nO " ^ i^n,m * Qm , (133) 

m=—oo 

where * stands for the usual convolution 

/ + CO 
de'fie-e')gie') , (m) 
-oo 

and Knm is the derivative of the logarithm of the effective S'-matrix between the strings of size n 
and m. In particular we have 

-Ko,o{e) ^ Ko{9) = ^^log^o'W , (135) 
for the interaction between physical rapidities, 

Zttz da & + i\n\/2 7r4c^^ + |n|^ 
for the interaction of 6 rapidities with Bethe strings of size |n| and 

2 

{9)= 2K2s{e) - Kn+m{e) + K\„i{e)6n^^ , n,m= 1,2,... , (137) 

for the interaction of two Bethe strings. Obviously i^n,m = if ?i x m < 0. 

It is interesting that even though these kernels appear as some quite complicated functions they 
all exhibit very simple fourier transforms Kn,m- More precisely we have 

and 

K„(u;) = e-HW2 ^ ^ = 1^2,... , (139) 
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so that the sum in (137) can be exphcitly done yielding 

Kn,m = coth (^M^ (^e~'-^\m-n\ _ _ ^^^^ ^ n,m = 1,2, .. . (140) 

A very useful formula for what follows concerns the inversion of the operator i^„,m when both 
indices n and m are restricted to be positive (or negative). In Fourier space 

S {5n,m+l + ^n,m.-l) ■, {n,m>0), (141) 

where the operator s (and its fourier transform) has the following form 

8(9) = |s(cu) = ^^l- (142) 

^ ^ 2cosh7r0 ' 2cosh(^) J ^ ' 

In particular we notice that 

Kq = 2s*Ki. (143) 

Having introduced all the necessary kernels we can proceed to construct the quantity of interest, 
the free energy at the temperature \/L. We have 

/(L) = min,„,^„ j dO (^Q^L cosh ttO - ^„ log (^1 + + log (l + ^ . (144) 

where we should minimize the integral by varying the densities of particles and holes keeping 
the Bethe equations satisfied. The first term in the integral is the energy density multiplied by the 
inverse temperature (which is L) and the second term given by the sum over n is the entropy density 
(see e.g. [17]). We use Bethe equations to write the variation Squ = —Sgn — Ylm=-oo ^'^nm * ^Qm so 
that the extremum condition 5f = yields a set of TBA equations 

oo 

= -en + Lcosh{7re)6nfl+ Yl -?^mn * log (l + e"^'") , (145) 



m=— oo 



where |2- = e ^" . The free energy evaluated at this extremum can then be written in terms of eo 
alone, so that 

Eo{L) = - Jy cosh(^0) log (l + e-^°(^)) . (146) 

The last two equations yield the finite size ground state energy of the SU{2) principal chiral field. 

Finally, we will show below that defining the incidence matrix Inm = Sn,m±i, the Y variables 
Ym = e^™- (m 7^ 0), Yq = e~''° and using the operator (141) these equations can be transformed 
into a local (in the discrete variable n) set of integral equations 

oo 

log {Yn) + L cosh(7r6') (5„,o = ^ InmS* log (1 + Ym) , -oo <n<oo (147) 

m=—oo 

and 

— cosh(7r0) log (1 + Yo) . (148) 
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It is remarkable that non-local kernels Knm disappear and at the end only a very simple kernel s 
appears in the final set of TBA equations. It is also remarkable, and still somewhat mysterious, that 
the SU{2)l and SU{2)ji wings are smoothly glued into one Y-system on a discrete set — oo < n < oo. 

To show (147) we should consider (145) separately for n > 0, n < and n = 0. Applying 
the operator (141) to (145) in the first two cases (it is convenient to rearrange them to have the 
combination 6mn + Kmn), we easily verify (147), except the case n = which we should consider 
separately. To find n = equation of the Y-system we consider (145) for n = — 1, 0, 1. The kernels 
Kn,m entering these three equations are 

-Ko,o = Ko = 2s*Ki , Ko,±i = Ki , K±m,o = -Km , m > (149) 

and, most importantly, 

iSm^l . (150) 

Thus if we convolute (145) for n = 1 with the inverse shift operator s and use that s * (Km+i + 
Km-i) = Km we get 

oo 

= -s*Ki* log (1 + e"'«) + ^Km* log (l + e"'™) - s * log (l 6+'^) . (151) 

m=l 

Notice that the last term is separated from the infinite sum because the m = 1 case in (150) behaves 
slightly differently than for the other m's. Moreover the sign of the exponent inside this log differs 
from that inside the logs in the infinite sum because we absorbed the first term in (145) into this 
last log. Similarly, for n = —1 we have 

-1 

= -s*Ki*log{l + e-'°) + i^i^i *log(l + e"''") -s*log(l-Fe+'-i) . (152) 

m=— oo 

These two equations can then be used to simplify the n = equation which reads 

= eo - L cosh(7r6l) + Kq * log (l + e^'°) - ^ K^m] * log (l + e"'™) . (153) 

m^O 

Indeed if we sum all these three equations we see that (i) the infinite sums completely cancel out, 
(ii) the convolutions with log (1 -|- e"*^") drop out as well by virtue of the identity Kq = 2s* Ki. We 
are thus left with the last terms in (151) and (152) plus the first two terms in (153) thus obtaining 
the last y-system equation (147) for n = 0. 

Moreover, for the functions g{9) analytic inside the physical strip lm(9) < 1/2 we have 

s*[g{6 + i/2)+g{9-i/2)]=g{e) (154) 

because 

<,(» + ,/2) + - ^ ^ r m ^ ^^^^ 



Ijl 2 cosh(7r(^ — x)) 2i Jpg sinh(7r(0 — x)) 

Therefore if Yn is non-zero inside the physical strip we can invert the s operator to find a set of 
functional equations, finally rendering the y-system for the PCF at a finite temperature 1/L 

Yn{e + i/2)Yn{e - i/2) = (1 + Yn-i{e)){i + y„+i(0)) (155) 

To fix a solution, this y-system ought to be supplemented by the large 9 asymptotics Yn — 
^-Sno Lcosh-ne i^j^g relativistic dispersion relation. Notice that these functional equations 

do not contain the dispersion relation explicitly: it appears only through the asymptotics of the 
y-functions. 
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B General solution in terms of Hirota functions 



In this appendix, we will give an alternative approach to the construction of solution for the 
energy of a general state of the SU (2) PCF in the periodic box. It will give a new NLIE defining 
the spectrum, different from the one of the section 3. This approach is the generalization of the 
approach we used in section 2.7 for the states with particles and non-trivial polarizations encoded 
in the Baxter polynomials Qu and Q^. 

As explained in the beginning of the section 3, for each solution to the Y-system equations, 
there are two natural solutions to Hirota equation which yield the same y's and are related by 
a gauge transformation. The expected analytic properties of these functions are described in this 
section. In particular we have 



rpV 


= g^r 


rr-iU 




1 


rpV 








= 9^9~ 





(156) 

which we will now use to completely solve our problem. First, as in the section 3, we find the gauge 
function g from the last relation. 



9 



le 



jLsmh(7ra;)/2 / ^ 



and plug it in the first two to find 



rpV 
^+1 
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-L cosh(7rx) 



-iLsmh(7rx)/2 / 



(157) 



-L cosh(7rx) 



($+$^$+$;r: 
($+$^$+$;r; 



(158) 



As in the section 3, we still have to relate T^'^ and but we do it here by a different relation. 
For that, let us define in the whole complex plane x the functions 



Fu{x) = (j){x 
and 

F^{x) = (j){x 



2-Ki 



4>{x) 
2-Ki 



T!^i{y)Quiy + i)/Quiy) 
{y-x + i/2)(t){y + i/2) 

TiMQv{y + i)IQv{y) 



dy + 



2TTi 



I 



dy + 



x) f T^MQu{y-i)/Qu{y) 

{y-x-i/2)(l){y-i/2) 

{y)Qv{y-i)/Qv{y) 



dy (159) 



rpV 
^+1 



dy, (160) 



{y-x + i/2)(l}{y + i/2) 27ri J [y - x - i/2)(l,{y - i/2) 

where the integrals essentially go along the real axis, but we should pass the contour in such a way 
that the zeroes Uj,Vj of Qu,v remain below the contour and the complex conjugated zeroes Uj of 
Qu,v stay above the contour. Using eq.(22) we Ccin show tlicit Tq ' cind. ^ are related to the values 
of the same analytic function v inside and outside of the analyticity strip, respectively: 



Fu{x) 



Imx > +1/2 
T^{x) \lmx\ < 1/2 , F^{x) 
Imx < -1/2 



Imx > +1/2 
T^{x) I Imx I < 1/2 



^^>+ Imx < -1/2 



(161) 
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Indeed, substituting from (22) T^^ = {T^"^ - $"Q")/Q"++ and the conjugate T^^ = {T^~ - 
^^Q^)/Q^ into the first and second terms of (159), respectively, we can separate the contribution 
of two integrals in (159) containing Tq into a single contour integral going around the physical strip 
(to realize it is useful to make shifts by ibz/2 for the integration variable). We shall also use the 
fact that ^latx— >iboo. If the point x is inside the physical strip we can contract the contour 
around the pole at y = rr and thus verify the middle relation in eq.(161). The other two integrals 
in (161), containing ^ and do not contribute since we can close the contours there around the 
upper and lower half-plane, respectively. The result is zero since there are no singularities inside 
by our previous assumption. The poles related to zeroes of Q^^v do not contribute since they are 
outside of these contours by definition. We should ensure by hand the analyticity of ^{x) close to 
the real axis 

J,U+QU ^ r^U^QU++ ^ ^ = , (162) 

which is the finite L deformation of the usual asymptotic axillary BAE for the n-roots. The relations 
for Fi, are found in a similar way. This equation can be shown to be equivalent to the eq.(70) derived 
in the main text. 

Notice that Tq and Tq are automatically analytic even slightly outside of the physical strip, 
because we can deform the contours to open further the physical strip. We will discuss this "extra" 
analyticity in the App.D. 

Finally, notice also that 

SO that (158) completely constrain the functions T"^ and out of which all other and can 
be written down using the resolvents Fu and F^. 

In the U{1) sector we have Qu = Qv = ^ and the two wings are obviously equivalent. We will 
have in this case = T^-^ and F„ = F^ will be given by (47). Equations (158) then reduce to the 
previously derived equation (51). 



B.l The main Bethe equations 

The main BAE reflecting the periodicity of the wave function and constraining the real zeroes 9j 
is given by (16): 

ri"(^,- ±f/2)r"i(0,- ±i/2) , , 

^ ' ' ' <^>u{0j±i/2)^ui0,±i/2) ^ ' 

Using the fact that T" has no zeros inside the physical strip and the denominator is regular for real 
6j we conclude that T^i{9j it i/2) ^ 0. This condition can be in fact interpreted as yet another 
form of the main BAE. It can be further simplified: using (156) we get 

~ ^ ^ ^ 



where we can use that, due to (22), for the n-wing 

iU -^U rj-\U-\- 



K = +T^^-T-t^. (165) 
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Substituting it into the eq.(164), and evaluating at 6j which is a zero of Tq we get 



l = y±!L_i| (166) 



or, using (22) for the v-wing, we get the simplest form of the main equation, easy to compare with 
the large L limit 

which was also derived in the main text eq.(67). 



C Proof of reality of T/ 



k 

In this appendix we shall analyze the reality of the T- functions. This is an important point to 
consider because the Hirota equation is solved explicitly by eq.(23) provided all are real. The 
goal of this appendix is to show that once the following equations (see (58)) are satisfied 

= S^^^ ±± (Igg) 

^« fQt+Qv- - QuQv 

= Ml^^ (169) 

fQt^Qu- - QvQu 

= g+g-<^>u, (170) 

and Tq* and Tq are real then all Tk are real and thus all the formulae in the main text go through and 
the y-system is indeed solved by (18). Before proving this statement we recall that eq.(66) follows 

rpU-\- rpU— rpV-\'rpV — 

directly from eq.(170) under certain analyticity assumptions and also from ^ — = '^^ ^ , 
which is a consequence of eq.(168) and eq.(169). Thus we can add the equation 

TS = 99T^ (171) 

to the equations at hand and proceed to the proof of the reality of the T- functions. 
Equation (168) implies 

- [-^- ' ^0 - ^/^i [— —) , (172) 

for some hi. Since Tq is real for real x the function hi is a real function. Eq.(169) implies 

l>„ = -i/i+(g"Q++Q-^ -5+Q,Q„) , Ti; = ih2{gQtQt -gQuQZ), (173) 

where again /12 is a real function. In virtue of (171) we have hi = h2 = h and by conjugating the 
first equation in (173) we find 

= ih- {9+Q--Qt+ -g-QyQu) = j-^a'g^^u , (174) 
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which means that the real function h{x) is periodic in the imaginary direction, h{x) = h{x + i). 
This in turn impHes that the function h" is also a real function because 

h{x - i/2) = h{x + i/2) = h{x - i/2) . 

Then it is simple to see that T" is real. We simply write, from eq.(172), 

r^u ^ + <^>u^ = ih- ( ^^^"^ - \ . (175) 

The reality of T^* is now manifest because the expression inside the brackets is purely imaginary 
and, as we have just shown, h~ is real. Proceeding in the same way one can see that all are 
real and thus the Hirota equation is satisfied by our solution. 



D Proof of analyticity of T_i in the physical strip 
D.l Analyticity in the U{1) sector 

From (51) together with i^(x it i/2 it iO) = -F(x it i/2 =p iO) — T-i{x) we have 

TiF++F~+ _ TiF+-F- 

^ [F++F — e-'^'^°''^('^^) + riF++ ~ [F++F — ^"^^ e^'^°^'^^'^^'> +TiF~- ' 

Using respectively the first /second equality we can smoothly move from real x into the upper /lower 
half complex plane provided Ti is analytic. In this way we can reach any x inside the enlarged strip 
|Im (x)| < 1 where Ti is regular. 

Notice that for large L we have 

( (t>^ ,+1/2 < Im (x) < +1 

T_i(x) ~ <^ ,-1/2 < Im (x) < +1/2 (177) 

( 0+ ,-1 < Im (x) < -1/2 

The denominator in the expression for T_i at x = Oj — i/2 is proportional to Bethe equations 
gi^iPjL _|_ ]^ — This is not a pole of T_i because the numerator at these points is proportional to 
To(%) = 0. 

However, for large volume, T„i could have poles at the analogue of the holes of the 9j BAE, 
close to the boundaries of the physical strip. 



D.2 General case 

In this subsection we will study the analyticity of T"]^(x) and Tf{x) for a general solution. We 
will show that for large enough L these functions are analytic inside the physical strip — 1/2 + e < 
Imx < 1/2 — € where e ^ when L oo. We start from eq.(158) and rewrite it as 

j^u ^ , ^-LcoshUx) ~ ^"iQm"^) ^uTli (178) 
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and similar for T^^^. Solving for T^-^ we get 

— L COsh{TTx) rpU— ^ 

Since is regular in the lower half plane and T^i is regular in the strip — 1 < Imx < 1 where 
($^^>~^+^~)*'^ is also regular the singularities of T^^ for —1 < Imx < could be only due to 
zeros of the denominator. 

As far as L or x are large for —1/2 < Imx < 0, the second term in the denominator is 
exponentially suppressed and to get a zero of the denominator we should be close to a zero of 
Qu- The points close to Uj where the denominator vanishes are in fact follows from 

the auxiliary finite volume BAE eq.(70). However these poles cancel with zeroes of the numerator 
rendering T_i regular at these points, a result we were familiar with already. Thus, we see that 
for large L the only poles that could appear must lie close to the border of the physical strip 
Imx = —1/2 where the exponent is oscillating. It oscillates faster for large x and we thus have 
poles condensing at infinity along the borders of the physical strip. 

We conclude that for the general solution - at least for large L's - the function T^^(x) is analytic 
inside almost the whole physical strip and could have poles only very close to the border. For small 
L it can probably happen that the singularities approach the real axe. That could indicate some 
singular behavior of the energy levels as functions of L such as the one observed in [5G] (see e.g. 
figure 10 in this work). It would be interesting to investigate these points in greater detail. 



E Details on Liischer formulae derivation 

If this section we shall present some details of the computation of the first finite volume correction 
to the asymptotic auxiliary Bethe equations, obtained by expanding (70) to the leading order (see 
section 3.2 for notation). We start by writing 

P{Uj) = exp [Ki{Uj - i/2 + iO) * pu - Ki{Uj + i/2 + iO) * pu] ■ (180) 

Notice that we introduced the iO's because P{x) was originally defined in (71) for x in the upper- 
half-plane. Removing the iO's by the use the Sokhatsky-Weierstrass formula we get 

PiU,) = exp [K^ * Pu - K+ * Pu + Yo/2] ^^^^ , (181) 

where the convolutions are understood in the principal part sense and Pu + Pu = log(l + Yq) — 
was used. Next we split the density pu into pu^ and pu^ as explained in section (3.2). The former 
contribution is purely imaginary and therefore it contributes to the exponent as 

A7 * pW - Kt * (-pI^)) = K2 * p^^^ . (182) 

Hence we finally obtain 

P{Uj) = exp (k2 * p(^) + (i^r * - c.c.) + ^ ) , (183) 
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with all convolutions understood as principal part integrals. It turns out that the first and the last 
terms in this exponent simply convert the Bethe roots Uj in (70) into their real parts, namely, 

^^ _ ^« + (« 



,exp i^2*P^^^(n.) = (184) 

riui)Qt+iui)^^^rYom\ ^ rm^t^m _ ^^^^^ 



The check of the first equality goes exactly as in (78) and we will therefore not consider it here. 
Let us explain how to check the second equality. Notice that this expression is equivalent to 



Next, we write the right hand side containing Yq ~ T^iTi/ (j)'^ as 



20+(^,)Q'jC/,) J \Qt+iU,)r{Uj)' • ^^^^^ 

The first factor in the r.h.s. is precisely Auj. This can be seen from expanding the second equation 
from (22) at x = Uj = Uj — 2Auj to leading order in Auj. Alternatively we can find the imaginary 
part of Uj by imposing regularity on the density (75) at x = uj. To simplify the second factor in 
(187) we write 



Q' 



(188) 



Evaluating the derivative of the numerator and using the leading order auxiliary Bethe equations 
Q--(?7,)<A+(C/,) + Qt^{U,)4>-{U,) ^ we find 

T^m^jUj) -a.logf^^l (189) 



thus identifying the second factor in the left hand side of (186) and completing our proof. Therefore 
the expansion of the auxiliary Bethe equation (70) simply reduces to (80), as announced in the main 
text. 

To simplify the Bethe equations (79) and (80) further we shall relate the convolutions in these 
expressions to particular derivatives of the 1"- function Yq. To compute these derivatives it is useful 
to notice that we can write Yq in terms of two simple pure phase functions and a„, 

Qw[x) 

as 

yo = K + lK)(a+ + l/a-) . (191) 

In this form, it is easy to compute the derivative of Yq with respect to Ok, Uk or Vk because we can 
use a simple identity 

de.ay,{x) = -TTiKQ{ei - x)au,{x) , dyo,aw{x) = 27riKi{wi - x)aw{x) . (192) 

Furthermore if we notice that the densities p^"^^ can also be simply expressed in terms of these new 
functions as 

pI'^ = at {at + , p1'^ = at {at + l/aZ) , (193) 

It is then a straightforward exercise to check the identities (83). 
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Figure 8: For several models the Y-system TBA equations for the ground state energy are known 
and can be represented by the diagrams such as the ones in this figure. Using the techniques 
developed in this paper it would be extremely interesting to compute their complete spectrum in a 
systematic way. 
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